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Photoelasticity and Aircraft Design 


H. T. JESSOP, T.D., M.Sc. and C. SNELL, M.Sc. 
(Reader in Photoelasticity, and Lecturer, Civil Engineering Department, 
University College, London) 


1. Introduction 


It is now more than 25 years since Coker and Filon first developed photo- 
elasticity as a practicable method of exploring the stress distribution in certain types 
of engineering components. All engineers must be familiar with the striking 
photographs of fringe patterns obtained by their method in stressed plates of various 
shapes viewed in a polariscope. In those early days observations were made on 
Xylonite models and the method was only applicable to components in which stress 
distribution was of a two-dimensional type. Moreover, with the apparatus then 
available investigations were somewhat difficult and many engineers at that time 
were very distrustful of the results obtained from what was, to them, a little 
understood laboratory technique. 


Since then, however, conditions have changed greatly in many respects. The 
discovery in 1938 of the “frozen stress” properties of some of the new synthetic 
resins opened the way to the application of the photoelastic method to three- 
dimensional stress distributions, the development of new materials has greatly 
increased the scope and accuracy of the method, experience has confirmed the 
reliability of results obtained, and engineers have become accustomed to the use 
of the techniques of the physics laboratory in their research problems. Today, in 
many branches of engineering, photoelasticity is recognised as one of the most 
valuable design tools, and new techniques are continually being developed for the 
application of the method to fresh problems. 


The Aircraft Industry is perhaps one of those in which the value of the method 
is least appreciated; this is probably because so many of the problems of aircraft 
design are of types which are unsuited to photoelastic investigation. For the 
photoelastic method, like all other practical methods of stress exploration, is not by 
any means of universal application. In some cases the shape of a component renders 
it unsuitable, while in others it may be found impossible to reproduce in a model 
the system of loading to which the prototype is subjected. Thus the thin panels 
which form the skin of wings or fuselage could not be investigated by this method 
because scale models of such parts made in photoelastic material could not support 
the necessary loads without buckling. It might be noted, however, that the trend 
towards thicker skins, if it should continue, would probably bring some of these 
parts within the scope of the photoelastic method. Again, a.considerable part of 
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the stress in some aeroplane parts is due to vibrations of the parts themselves, and 
it would be impossible to freeze such dynamic stresses into a photoelastic model. 
Electric resistance strain gauges applied to the actual structure are very properly 
used for the exploration of stress in such parts. 


There are, however, many aeroplane components eminently suitable for photo- 
elastic investigation; components whose form is compact and for which the type 
of loading applied in service may be very clearly defined. Examples of such 
components may be found in the short spars and spar attachments connecting 
sections of the aircraft, in mounting brackets for engine or undercarriage assemblies, 
and in hinge brackets taking the reactions from control surfaces. All of these are 
vital to the safety of the aircraft, all are from time to time subject to failure, and in 
none of them is the stress distribution at all well known. These are, moreover, parts 
which are not well adapted to exploration by strain gauges, for the stress 
concentrations in them almost invariably occur in fillets and around holes and 
cut-outs where the use of gauges is extremely difficult. Granted that the loads to 
which such parts are subjected vary in both magnitude and direction, and that the 
failures are almost always due to fatigue—neither of which conditions can be 
reproduced in a frozen-stress model—but fatigue failures are very closely related to 
the peak values of the stresses at points of greatest stress concentration, and these 
points may be identified, and the stresses for any specified load distribution assessed, 
by the examination of a photoelastic model under a system of static loading. The 
photoelastic method therefore provides a rapid and simple means of comparing the 
relative strengths of different designs in such components for any defined load. 


2. Outline of the Method 


A scale model of the part to be examined is first made in one of the optically 
active plastics. This model may be machined from a solid block or cast in a mould, 
the method chosen depending upon the plastic employed and the form of the part. 
The model is then subjected to a system of dead weight loads reproducing 
qualitatively one of the systems of loading known to occur in service. It is then 
placed in an oven and heated to its “critical temperature ” (varying from 85°C. to 
120°C. for different plastics) and allowed to cool slowly. After cooling the loads 
are removed and the stresses are then found to be “frozen” into the model. 


One method of investigating these models is to cut thin slices from them and to 
examine the slices in the polariscope using the two-dimensional technique. In 
order to localise the observations these slices are cut as thin as possible (a thickness 
as small as 0.04 in. can be obtained without the machining process affecting the 
frozen stresses), and the optical effect upon light passing through a slice is not large. 
The close fringe patterns, therefore, which generally characterise the two-dimensional 
problems are not obtained, the maximum effect at a point of stress concentration 
rarely producing more than about the first “fringe,” but modern apparatus and 
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methods allow measurement of even this small effect to an accuracy generally better 
than one per cent. Examination of such a slice will yield the mean value through 
its thickness of the difference of what are called the “ secondary principal stresses ” 
in the plane of the slice, and if the slice is sufficiently thin such a mean value will 
give a good approximation to the actual value at any point. The determination of 
the principal stresses from observations of these secondary stress differences, 
however, can only be achieved in practice at points where the direction of one 
principal stress is known, and it is fortunate, from the point of view of practical 
application of the method, that the stress concentrations in which the designer is 
interested almost invariably occur in regions in which this is the case; that is, on 
the free surface of the component or in a plane of symmetry. 


A second method of examination, the technique for which has only recently been 
developed, is known as the “scattered light” method. This consists in passing a 
very narrow beam of polarised light through the solid model and viewing it from a 
direction perpendicular to that of the beam. For this method the model must be 
made in a photoelastic material which has the property of scattering light well. This 
is ensured by incorporating in the material during casting a small quantity of very 
finely divided solid. The material then has a slightly misty appearance, and the 
path of the light through it is clearly visible as a narrow illuminated pencil. When 
there are stresses in the model this pencil is seen to be crossed by dark interference 
bands or fringes, and if a graph is plotted of the positions of these fringes, the slope 
of the graph at any point gives a measure of the optical effect of the stresses in the 
material at that point. Thus in this method, instead of measuring the mean effect 
of the stresses through a slice of finite thickness, a measure of their effect at any 
point on the beam is obtained and, since the beam itself may be of a width less than 
0.01 in., the method gives a high degree of localisation to the observations and may 
be used with quite small models. 

This method is subject to the same limitations as the slicing method in that the 
stresses recorded at any point are only principal stresses when the direction of the 
light beam is that of one principal stress at the point. This condition, however, can 
always be fulfilled in the most important region, i.e. at a point on the free boundary 
of a model, by allowing the beam to enter the model normal to the surface. 
Measurement of the fringes seen in such a beam enables a graph to be drawn, and 
the slope of the graph at the point corresponding to the point of entry of the beam 
gives a measure of the difference of the principal stresses in the plane of the surface 
of the model at that point. The approximate directions of these principal stresses 
are also determined by observing the inclination of the model necessary to give the 
sharpest definition of the fringes, so that when the slicing technique is used at a 
later stage to measure the separate stresses, the direction in which slices should be 
cut is known. Observations by the scattered light method are also used for 
preliminary exploration of the surface stresses in cases where the location of the 
points of stress concentration is not immediately obvious. This enables the 
investigation to be made without having to cut a large number of slices. 
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There are other methods of examination which may be needed in certain special 
cases, but for most components of the type mentioned earlier, a combination of the 
two techniques which have been outlined here will give all the information the 
designer needs. 


3. Two-Dimensional Analogies 


While the accurate determination of the stresses in a component may require 
the use of the three-dimensional techniques outlined, there are many cases in which 
a large amount of useful qualitative information may be obtained from the 
examination of models which reproduce in two dimensions some of the characteristics 
of the actual component. Such models may be cut quickly and inexpensively from 
polished sheets of suitable photoelastic material, subjected to load, and examined 
directly in the polariscope. This preliminary investigation is frequently useful in 
enabling a qualitative comparison to be made between various modifications of a 
design, for although the stress distributions in such plane models are not identical 
with those in the prototypes, they indicate with a considerable degree of reliability 
which type of modification is likely to be most effective in reducing dangerous 
stress concentrations. 


In some cases it is possible that these two-dimensional analogies may yield 
evidence which is sufficiently conclusive to make the examination of a three- 
dimensional model unnecessary. In others it may not be possible in the model 
to reproduce conditions sufficiently close to those in the prototype, and a fuller 
investigation will then be needed. 


4. Illustration of the Photoelastic Method 


As an example, an investigation has been made of the stress concentrations in 
an elevator hinge-bracket from a Spitfire. While this component may not be of 
special current interest, as it is understood that this particular design is now obsolete, 
the stress problems involved in its design are of the same type as those occurring in 
similar components of more modern aircraft, and the methods described here would 
be equally applicable to such components. 


The loading of this particular bracket in service is very simple and is indicated in 
Fig. 1. The load (estimated at 600 Ib.), is applied through a steel pin fitting in a 
bushed hole in the lug of the bracket, and acts in the plane of the lug. The base of 
the bracket is secured to the channel-section rear spar of the tailplane by four 2 B.A. 
bolts, recesses (A) 4 in. in diameter, being machined in the base to provide 
parallel bearing surfaces for the securing nuts. The bush (B) is secured to the lug 
by four countersunk rivets through its flange. 
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A 600 Ib. 


La 


Slice at Xx’ 


Fig. 1. 
Drawings of bracket. 


Inspection of the design suggested that the two regions of greatest stress 
concentration would be found, (i) at some point on the boundary of the hole where 
the load is transmitted from the bush to the lug, and (ii) in the neighbourhood of 
the two securing bolts which are in tension under the load. There were indications 
that the stress in the second region would prove to be the more serious, for the 
sharp re-entrant angle produced by machining the recesses would cause a further 
concentration of tension very close to the area of load application under the nut. 
A preliminary two-dimensional test was made to study the effect of variation of 
design in this region. 

Figure 2(a) shows the fringe pattern in a plate cut to represent a section of the 
bracket through the two upper bolt-holes. The model was lightly secured to a fixed 
horizontal plate by two clamps, the jaws of which correspond to the nuts in the 
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machined recesses, and a vertical load of 10 Ib. was applied to the vertical arm. 
Fig. 2(b) represents in two dimensions the effect of providing bosses which form 
seatings for the nuts and which may be machined without introducing a re-entrant 
angle. The load on this model was 12 lb. An inspection of the fringe pattern shows 
that while the 10 Ib. load in the first model produced more than three fringes in the 
angle, the 12 Ib. load in the second model produced less than three fringes. A more 
accurate measurement gave the ratio of the peak stresses per unit load in the two 
cases as 5 to 3. 


The fringe pattern in Fig. 2(a) is characteristic of a sharp angle under tension. 
At first sight it would appear that the boundary stress is higher on either side of the 
angle than it is at the vertex, but examination of a much enlarged image shows that 
the highest stress actually occurs at the vertex, but is confined in that region to an 
extremely narrow strip along the boundary (see Fig. 7). 


For the three-dimensional investigation a full-size model was cast in “ Marco” 
resin using a “ Vinamold” mould, with an actual bracket, supplied by Vickers- 
Armstrongs, as pattern. The bush was turned from another block of resin and 
fitted into the lug after casting. Fig. 3 is a photograph of the finished model in its 
loading frame. The inverted loading was adopted to allow for the model being 
immersed in a tank of glycerine and water during the “freezing” process, this 
precaution being necessary to prevent the formation of a “rind” stress in the 
material due to loss of moisture from the surface on heating. 


Special care has to be taken in the method of applying loads when freezing 
stresses into these models, especially when it is required to investigate the stresses 
near to load points, for the yield of the photoelastic material is large compared with 
that of the material of the prototype, and this may lead to serious changes in load 
distribution if suitable precautions are not taken. In this model the method of 
bolting the bracket to its support needed special attention. To ensure that local 
yield of the plastic round the bottom bolts did not allow the base of the bracket to 
lose contact with its support, these bolts were provided with spring washers 
(P, Fig. 3) whose strength was just sufficient to maintain contact under the load. 
Also, thick washers of the “ Marco” resin were inserted under all the nuts to prevent 
the edges of the nuts cutting into the softer material of the model. One feature of 
the prototype which was not reproduced in the model was the riveting of the bush 
flange to the lug. To obtain a satisfactory reproduction of this, using rivets of the 
same material as the lug, would have necessitated making a larger scale model, and 
since the observations on the first model showed conclusively that this was not the 
critically stressed region, it was thought unnecessary to undertake this. 


It may be noted, however, that the way in which the load is shared between 
rivets and bush in such a construction as this might be important in other similar 
components. The distribution will obviously be related to the amount of initial 
tension in the rivets, and tests would have to be made using some method of varying 
this factor. This problem might well be the subject of a future investigation. 
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Fig. 2. 
Fringe patterns in two-dimensional models. 


Fig. 8. 
Enlarged photograph of fringes in corner of slice cut along XX’ (Fig. 1). 
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Fig. 3 
Photoelastic model in its loading rig. 


Fringe pattern in lug. 
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Photoelastic model in its loading rig. 
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The model was subjected to a load of 3 Ib. and heated in its glycerine and water 
bath to 90°C., at which temperature it remained for half an hour before being 
allowed to cool to room temperature in about 15 hours. The load was then removed 
and the model taken from its frame for examination. The bush was first removed 
and the lug of the bracket examined in the standard polariscope. Under the loading 
conditions of this test the greater part of the lug has approximately a simple two- 
dimensional stress distribution, and the principal stress differences may be deduced 
directly from the fringe pattern seen through the thickness of the lug. Fig. 4 is a 
photograph of these fringes taken in monochromatic mercury green light, and Fig. 5 
shows the distribution of lines of principal stress in this region as deduced from 
the isoclinic fringes. 

At points on the “free” boundary the fringes indicate the actual magnitude of 
the boundary stress (the stress normal to the boundary being zero), and this is seen 
to be a maximum at A (Fig. 4) where the stress in the model produces 1.3 fringes 
through the thickness of the lug. The relation between optical effect and stress is found 
from a calibration test-piece cast from the same mix as the model and heated with 
it under a simple load. From this the “fringe value” of the material is found. 
This is the stress difference in Ib./in.? which would produce one fringe in light 
passing through one inch thickness of the stressed material. For the material of this 
model the fringe value was found to be 4.47 Ib./in.* per fringe per inch, and from 
this the tension at A was calculated as 19.3 lb./in.* in the model, which corresponds 
to approximately 3,900 Ib./in.* in the prototype under a load of 600 Ib. 


On the boundary of the hole at B the stress difference indicated by the fringe 
pattern is the algebraic difference between an approximately tangential tensile stress 
and the compressive stress due to the load transmitted by the bush. The 2.86 
fringes recorded at B correspond to a stress difference of 9,600 Ib./in.? in the 
prototype, or 4,800 Ib./in.? shear stress. To obtain the magnitudes of the separate 
stresses in this region, observations by the scattered light method were made. The 
end of the lug was cut off to allow a beam of polarised light to be passed into the 
model normal to the boundary of the hole at various points around the portion AB 
of the circumference (Fig. 6). The graph of the fringes at the point C where the 
circumferential stress was found to be a maximum is shown in Fig. 7. The slope 
of the graph at O shows a fringe gradient of 6.5 fringes per inch, and this 
corresponds to a tension in the boundary of 29 Ib./in.? in the model or 5,800 Ib. /in.? 
in the prototype. The normal thrust due to the bush is then deduced as 
3,800 Ib. /in.? 

To examine the stresses in the region of the securing bolts a thin slice was cut 
from the frozen model, the mid-plane of the slice corresponding to the section across 
XX’ (Fig. 1). The slice was 1/16 in. thick, and when examined in the ordinary 
polariscope showed less than one fringe over most of its area, although at least one 
fringe was discernible in the sharp angle corresponding to the inner edge of the 
recess A. When examined with a magnification of 15 times, however, the second 
fringe was seen at the extreme point of the angle. Fig. 8 is a photograph taken with 
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600 Ib. 


Fig. 5. 
Lines of principal stress in lug. 
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Path of 
Light Beam. 


Compensator 


Fringes seen in 
Light Beam. 


Fig. 6. 
Method of examination by scattered light. 
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Fig. 7. 
Graph of fringes in scattered light. 
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this magnification. Thus the stress concentration at this point produces two fringes 
through a thickness of 1/16 in. and this corresponds to a stress of over 28,000 
Ib./in.? in the prototype under 600 Ib. load. Examination of the slice in a tilting 
Stage microscope confirmed that the plane of the slice is a principal plane in the 
region of the angle, and this stress therefore represents the true maximum tension 
in the corner. 


The original bracket was then modified by the addition of four cylindrical plugs 
which fitted into the machined recesses and projected slightly to form bosses on which 
bearing surfaces for the nuts could be machined. A new model was cast from this 
and stresses were frozen into it as before. The modification produced no change 
in the stresses around the bush, but examination of a thin slice corresponding to 
that cut from the previous model showed no trace of the high stress concentration. 
The optical effect through a thickness of 1/16 in. was too small to produce the 
first fringe, and so gave nothing which could be photographed. Measurement by a 
compensation method showed that the maximum stress in this region produced about 
0.9 fringe, corresponding to a tension of about 12,500 Ib./in.? in the prototype. 


5. Conclusion 


The investigation described, while by no means comprehensive, has shown the 
capabilities of the photoelastic method in obtaining information as to the stress 
distribution in such components as this. In this example it shows that while the 
lug of the bracket is amply strong, there is a very high, although localised, stress 
concentration in the sharp corners of the machined recesses in the base. How 
the material of the bracket will react under this concentration is a question for the 
metallographist, but it seems likely that repeated application of the specified load 
would lead to the formation of fatigue cracks and ultimately to failure. In any case 
the existence of such high concentrations is obviously undesirable. It may be noted 
also that the abnormally high sensitivity to fatigue shown by some of the new light 
alloys makes it more than ever necessary for the aircraft designer to know the 
magnitudes of the stress concentrations in all components. The photoelastic 
method is at present the only means of determining exactly how high such 
local stresses are. 
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A Note on the Design of Ducted Fans 


BRYAN THWAITES, M.A., Ph.D., A.F.R.Ac.S. 


Summary: The design of ducted fans, and the estimation of their performance 
for conditions other than those of the design, are usually based on the equations of 
the so-called strip-theory. Existing methods of solution invariably employ some 
form of iteration which is more or less rapidly convergent according to the experience 
of the computer and the detail of the method. 


In this note, a solution is given which, although it is strictly an iteration, is 
explicit in the sense that in practice no iteration is necessary. It is shown how both 
design for a given performance and performance for a given fan can be calculated 
in one direct step. 

Conditions for maximum efficiency are considered in an. Appendix. Here a 
suggestion of R. Baldur’s is used to give a simple approximate solution of the 
strip-theory equations which can be developed to give conditions for maximum 
efficiency, for the special case of a constant pressure rise through the fan. 


Notation 
Suffix “0” denotes an approximate value (except for r,, a, and z,). 
r radial distance 


r, root radius 


R tip radius 

q 

duct 

Q torque 

T thrust 

P power 

© angular velocity 

N number of blades 

c chord 

F flux through fan per second 


Paper received November 1950. 
[The Aeronautical Quarterly, Vol. III, November 1951] 
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p air density 


axial velocity 
Pp pressure difference 
a, slope of two-dimensional lift curve 
@ incidence 
—a,  no-lift angle 
lift coefficient 
Cy drag coefficient 
y  tan-?(Cp/C,) 
a, inflow factor 
n efficiency 
n, efficiency of blade element 
6 blade angle 
6-2 
8 defined by equation (12) 
8 defined by equation (19) 
G,H defined in equation (1) 


1. Strip Theory Equations 


Let P, T, N, © and » be respectively the motor power, thrust, number of blades, 
angular velocity and efficiency of a fan which occupies the annular space r, =r<R. 
Let u, c, a, 9, 2, p, Cy and Cy be respectively the axial velocity, chord, axial inflow 
factor, blade angle, incidence, pressure difference, lift coefficient and drag coefficient, 
these eight quantities being functions of r. Then if (dT/dr) ér and (dQ/dr) dr are 
respectively the thrust and torque due to the annulus between r and r+4r, the strip 


theory equations may be written 


=4p New? Hrcot¢ 
sing cos? sing cos¢ cosy’ sing cos? sin¢@cos¢cosy 


dr’ rQ(i-a,) 
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A NOTE ON THE DESIGN OF DUCTED FANS 


Before showing the method of solution of these equations to give certain required 
results, it may be shown briefly how the requirements of the fan. determine the values 
of some of the parameters at the beginning. 


There are two usual methods of approach. In the first, the power of the fan 
and the flux through it are known and it is desired to design a fan of the maximum 
efficiency. This involves determining the values of u, R, r,, p and © to give 
maximum efficiency, which is an impossible task generally, although it may be done 
toa good approximation in the special case p=constant, as is shown in the Appendix. 


The second approach is the more usual in which some of u, R, r,, p and © are 
fixed at the beginning. In the choice of these values, previous experience will 
naturally be used to avoid values which might lead to a particularly low efficiency. 
There are many possible procedures but the following are points which usually arise: 


(i) u is assumed known and also the required down-stream values of pressure 
and velocity. This implies a knowledge of the power output of the fan and hence, 
by using a reasonable value of the efficiency, of the power input. This leads to a 
choice of motor of which the power P is fixed (but the angular velocity will be 
determined from the following criteria). P is therefore taken as a fixed quantity. 


(ii) R will probably be fixed by non-aerodynamic considerations. The aero- 
dynamic conditions on r, are that the C, at r, is allowable, with possibly a condition 
for non-interference between the blades and that r,/R should not be too small for 
the sake of efficiency. 


(iii) An upper limit on {2 is fixed by the Mach number at the tips, r=R, and 
by the centrifugal stresses which may also impose condition on c. 


(iv) N and ¢ are fixed in conjunction to ensure allowable values of C; and 
sufficiently high values of Reynolds number, while ensuring non-interference. 


More detailed information concerning these preliminaries may be found in 
previous papers on fan design. For the present purposes, the above considerations 
enable us to regard P, ©, N, u, as fixed, and a good approximation (according to 
the accuracy of the efficiency estimate) of p known. 


The value of 4, the blade angle, is required. 


2. General Process when u and p are given arbitrarily 
From equations (1), 


4T _ pNeu*. cose+y) 
 Qnrdr —“sin® cosy 
(2) 


and tang= ro {1+ sing COS¢ Cosy 
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Elimination of C;, between these two equations gives 


rQ 
— =coto+ 
<p 


tan(9+/y), 


which can be rearranged as 


{8 — [cos B {4 cos y( 1+ +siny } 


u p 


The desired form of p but not its scale, will have been determined. Now 


R 
and on a reasonable value, »,, of » being chosen, a first approximation to p, p=p,, 
may be made which satisfies (6). The error in p, will be one or two per cent. and 
it is easily verified that in usual cases the change in ¢ due to that possible change in 
p, from (4), is practically insignificant. 
Similarly it can be verified that small changes in y do not affect the value of 
given in (4). y is regarded as a function of C, (and Reynolds number) and an 
approximation C,, to C, is obtained from equations (2) in the form 


Cro Po sin %o tan = 


= 9) (7) 


All quantities on the right hand side of equations (4) and (5) have now been 
found to a degree of accuracy sufficient to determine the final value of 9. 
Before going on to calculate C;, and so 2 and 6, it must be ensured that the 
torque is Q. If suffix “0” denotes an approximate value, 
R 
P,=2.Q,=270 p,r*tan(o+y)dr (8) 
"0 
and the better the choice of 7,, the closer will P, be to P. If now p is taken equal 


to (P/P,) p,. ¢ will be unaltered, the torque will be exactly Q, the efficiency will be 
n=n, P/P, and the thrust will be exactly »P. Thus C, is given from (2) as 


_ 4rrp, P sin’ ¢cosy (9) 
~ pNeuw P, cos(¢+y7) 


Ci 


from which and finally 6=9+ follow. 
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The design has been accomplished therefore in a direct manner. The only 
awkwardness is the inevitable numerical integrations for Q, and P, and one of these 
can be avoided in the special case u=constant, p=constant, considered in Section 3. 
An arbitrary initial choice of 7 is also inevitable, but a good choice can usually be 
made and it is noteworthy that the error in the initial choice does not involve any 
recalculation. 


3. The Usual Case, u=constant, p=constant 
In this case, the work can be considerably simplified and a direct estimate of 
», made at the beginning. 
nP 


where the area of the duct is A=x(R?—r,”). Substitution for p in (4) and (5) gives 


o= { -y+sin-*| cos 8 { “(14 eos y+siny} 


] 


y is given sufficiently accurately from the approximate value of C,, given, from (7), by 


4x nP sin 9, (13) 


Ci. = 
pNe® Au? rQ 


» may be estimated as follows. From (1), and (10) dQ/dr=[27»Pr/(uA)|r tan (¢ +7). 
Now rtan(@+y) does not vary greatly across the duct, and if 


rtan(@+y)=L,+L, 


where 


then Q= [1+O(R-r,y] 


in which the implied constant in the O is itself small. A good approximation to the 
power P is therefore (yPQ/u) [rtan(¢+y)]--,, and the approximation, »,, to 7 is 
therefore chosen as 


in which a guessed value of » may be taken in the calculation of ¢. 
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¢ can now be calculated from (11), and then Q, is given by 


2n "o 2 
r tan(¢+y)dr. 


The true efficiency is then »=», P/(QQ,) and so C, is given by (2) and (10) in 
the form 


_ _P__ sin* ¢cosy 
pNcAuw® 2Q, cos(¢+y)* 


(15) 


From this, z and lastly 6 can be evaluated. 


4. A Worked Example 


A ducted fan is required to absorb 15 h.p. in a duct of area 11.5 ft.? in which 
u=77 ft./sec. The motor to be used runs at 1,500 r.p.m., r, is chosen reasonably 
as 4R, and this fixes R. 27r,=7.47 ft. and so Nc=3 ft. is taken. The tip speed 
is 353 ft./sec. and so there will be no compressibility effects. From (13) it is found 
that C;, varies approximately from 0.55 at the root to 0.22 at the tips. An aerofoil 
section with a design C, approximately equal to 0.4 and a low-drag range of, say, 
0<C, <0.8 is suitable, and N.A.C.A. 65,-418 is chosen. To obtain a low drag 
coefficient, the section Reynolds number should be at least 10°, and so N=3, c=1 ft. 
are taken. The efficiency », is given by [u/(r, 2)] [cot (¢+-y)],_, ; 7, =1.72 ft., and 
with taken as 0.87 in this calculation for ¢, and y as 1.30 to correspond with a 
C, of about 0.4, n, equals 0.89. C,, is then given from (13) and using the value of 
the Reynolds number, y, can be written down. With »=n,=0.89, 8 and ¢ follow 
from (12) and (11), and then dQ,/dr. Integration of this gives the torque Q, and 
so the power P,=0Q, which, it is noticed, is near 8,250. Cy, is then calculated 
from (15) and then 2 and @. It is not necessary to do any more calculation since the 
whole process is self-checking. The input power is now 8,250 ft. Ib./sec., and the 
efficiency is », P/(Q Q,)=0.892. The whole work is set out in Table I. 


5. The Inverse Problem 


A fan having been designed for one condition may operate at other conditions 
and it is desirable to be able to calculate speedily the case of 6 given. A relation 
between Q and {2 must be assumed (the motor characteristics for a fixed 
control setting). 
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TABLE I 
p=8,250, A=115,  u=77, 0=157, R=225,  r,=1.189 
,=1.719, N=3, e=10, <—~=03904, 

2pAu' 
) in R.N. dQ, 
Is) 1189 0.4125 127 046 1.04 16.24 2360 33.78 0517 209 257 


1.366 0.3590 143 O41 1.13 1422 2054 3862 0446 1.39 21.9 
1.543 0.3179 160 0.37) 1.23 12.65 18.19 43.72 0.393 0.85 19.0 
1.719 0.2853 1.76 0.33 1.33 11.39 16.34 49.05 0.352 045 168 
1.896 0.2587 1.93 0.30 143 10.35 1482 5460 0.319 0.12 149 
2.073 0.2366 2.10 0.28 155 948 13.55 6041 0.291 -0.16 13.4 


ich 2.250 0.2180 2.27 0.26 1.63 8.75 12.49 66.33 0.268 —040 12.1 
bly 
sed Atr=r,, if »=0.87, y=1.3° Q,=52.45, 02 Q,.=P,=8235; 
ind B=11.50°, ¢=16.33°, n, =0.89. P/P,=1.002; n= 0.892. 
ay, 
rag From equations (1) 
ft. 
ind 
= tan (14 and ¢=6—2 . (16) 
sin COS @ 
of 
nd 
“ is assumed, —a, being the no-lift angle of the blades, elimination of H, 2 and C, 
* from (16) and (17) gives, after some working, 
+4sin~' { cosa{ 1 + —1+ cos (6 + z +») |} (18) 
8zrcos y 

where tand=cot(#+ 2, —y)-4 (14 ) Nea. cosec (6 + 2, — y) (9) 
ns 
on Equation (18) therefore gives ¢, and a and C, follow immediately. dQ/dr can then 
ed be evaluated, and hence Q. This calculation is performed for a few values of °, 

so that a value of 2 can be chosen for the fan’s performance to suit the motor. 
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APPENDIX 


Equation (4) gives ¢ in terms of u, r, ©, p, p and y. 


The author is indebted to R. Baldur for pointing out that a rather simpler 
formula can be obtained for ¢ by using the local blade efficiency, »,, as a parameter 
in place of y. For from equations (1), 


H 1dQ/daT_ u 


tan 2 
+Y G rdr dr rQ (20) 
and so equation (3) can be written as 
rQ [ p ] 
ote= 1- | - 2 
u 2r? QO? pn, (21) 


For this formula to give as accurate a value of ¢ as equation (4), a good estimate 
of 7, must be made. This value of », inevitably involves y, which, however, it is 
easy to estimate. Elimination of between (20) and (21) gives 


() 
tan y( - =) te 


n,=1- 
rQ ou 


OF p 
which is true to O (y?) and O[p?/(4r* 04 p?)].. and p/(2r? 0?p)* rarely take values 
greater than 0.1 and are usually much less, and so approximations to the order of the 
square and product of these quantities are likely to be good enough. Thus, using 
(22), (24) and (20), a solution to this order of approximation can easily be found 
of the basic equations (1); in particular, equations (8) and (9) become 


R 

200. | — +o ] 

Q,=2 r| 14 OF p +tany dr (23 

4x p,u? P ‘+ ao wid 
G. Nc P, 2r? 0? p | (u? +7? 2?) ( 
0? 


(24) therefore gives C; and so z and 6=#+a directly from p,. 


The original formulae (4) to (9) do not involve much more computation than 
those of this Appendix and are probably preferable by being exact. However, the 
approximation given here enables the integration in the case p=constant, 
u=constant, to be performed, a direct estimate of the efficiency to be made and 
conditions for maximum efficiency to be formulated. 


*n, =[1/(2J)] [J+ + Oy), where J=POp/(2p). Therefore, J>1 and so must 
be less than 4. 
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If p is constant, and a mean value for y is taken, 


= 1 2yu p R 
whence n= 1- (R 37, OR )+ log (25) 


In any problem the power P will be fixed, and also probably the flux 
F=7u(R?—r,*) of air through the fan. (25) can therefore be rewritten as 


1) 


q=r,/R 


and the problem now is: given F and P, what values of R, g, y and give a maximum 
value for »? First, the smaller y, the greater 7. A clue to the answer is given by 
taking y=0, when R©) must be as big as possible (and the limitation will be imposed 
by the tip Mach number) and g as near unity as possible. Thus if g=1—A and 
working to O(A*) in », the value of g which makes » a maximum is 


F (26) 


f2 PR= 
q=1- ‘ 


and the resulting expression » shows that both R and {2 should be as large as possible. 
It should be remembered that y has taken a constant value, but this is consistent with 
the analysis being to O(y*). For high efficiency, the design consists therefore of 
choosing large values of R and {, and of then finding p and hence r, from (26). An 
incompatibility may arise if u is also fixed by the design, for F, gq and u are related 
by F=xu R?(1—q?): if the value of q given in (26) is incompatible with this, a 
value as near as possible should be chosen. 
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A New Treatment of the Tuned Absorber 
with Damping 


D. C. JOHNSON, M.A. 
(Department of Engineering, Cambridge University) 


Summary: The design of the damped tuned vibration absorber, which has 
previously been treated algebraically, is studied by the use of a three-dimensional 
admittance diagram. The advantages of this method are that the basic principles can 
be more easily understood and that it is not necessary to simplify the vibrating 
system so much. The results obtained by the use of the method show that under 
some conditions the algebraic treatment does not give the optimum tuning, but 
that the differences in performance of dampers which have been designed by the 
two methods will not be large. It is shown how the method may be extended to 
more complex problems. 


1. Introduction 


The damped tuned absorber is a device which is used for reducing torsional 
vibration in engine crankshafts. It is convenient to discuss the problem of damper 
design in terms of the idealised dynamic system which is used for calculating the 
natural frequencies of torsional vibration of crankshafts. Such a system is shown 


I 


| 
I 
| 
| 
| 
Fig. 1. 


in Fig. 1 to the right of the broken line. It consists of a number of rigid discs, which 
have moments of inertia J,, J,, .... and so on, and which are connected by light 
shafts having torsional stiffness A,, A, .... and so on. The discs will be acted on by 


Paper received February 1951. 
{The Aeronautical Quarterly, Vol. III, November 1951] 
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torques P, e't+*), P,e'+*), .... and so on, respectively, due to the gas and 
inertia forces on the pistons. These torques produce torsional vibration of the 
system. The absorber consists of an additional disc, moment of inertia Jp, attached 
to the crankshaft by a member of torsional rigidity Ap; some arrangement for 
dissipating energy is connected across this member so that energy is absorbed if there 
is relative motion between Jp and 7,. Such an arrangement may use dry or fluid 
friction, or it may consist of a rubber cylinder which absorbs energy because of its 
elastic hysteresis. The absorber is represented to the left of the broken line in 
Fig. 1. In this paper it is assumed that the energy absorbing mechanism has a 
linear characteristic, so that the torque T> is related to the relative angular velocity 
of J, to 7, by an equation 


Tp=1 
where u is a constant, and 4, is the displacement of Jp. The problem is to choose 
values for Ip, Ap and » which will prevent excessive torsional vibration of the 


crankshaft throughout the range of variation of , the pulsatance of the 
applied torques. 


Notation 
The numbers and suffixes 1, 2, .... refer to the discs. 

I moment of inertia 
torsional stiffness 

P amplitude of applied torque 

®  pulsatance of applied torque 
@ phase angle of applied torque 

Ty absorber torque 
damping constant 

a,, direct admittance at 1 

@,, cross admittance from 1 to 2 

6, amplitude of motion at 1 

@ absorber admittance 

X and Y _ real and imaginary parts of 1/2» 
quality factor for absorber 


Q 
K ratio of applied frequency to damper’s natural frequency 
F force exerted by the system on the absorber 

P, 


force applied at 2 (in Section 9). 
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2. Previous Treatment of the Problem 


An algebraic solution of the problem is given by Timoshenko"). In order to 
simplify the mathematics he replaces the dynamic system of Fig. 1 by that of 
Fig. 2(@). He then shows that, if values of Jp and Ap are chosen, a family of curves 
relating amplitude of J with may be drawn for different values of u. All the curves 
of the family pass through two points; that is to say, there are two values of © for 
which the amplitude of J is independent of the damping. The surprising result may be 
used for choosing the values for the constants. Jp and Ap are chosen so that the 
two particular amplitudes are equal and » is then chosen so that the amplitude. 
frequency curve has these values as maxima. 


I, I I 
! 
i, 
| 
on 
7 
Fig. 2(a). Fig. 2(5). 


It is stated by Zdanowich and Moyal) that Timoshenko’s treatment is not 
always satisfactory, presumably because of the simplification involved. This 
simplification is not only the replacement of the dynamic system by one which has 
a single degree of freedom; it also involves the assumption that the torque is applied 
at the point where the absorber is attached and that the amplitude at this point is 
the criterion of the severity of torsional vibration. The design method used by 
Zdanowich and Moyal does not require any of these assumptions to be made. The 
reason for the development of the treatment given in the present paper is that a 
method was required which would help engineering students to understand the 
fundamental ideas upon which any full treament, such as that of Zdanowich and 
Moyal, must be based. 


3. Outline of New Treatment 


The treatment depends upon the use of mechanical admittances (Duncan”’). 
The principle is to replace the algebraic problem by a geometrical one; this is done 
by considering certain characteristic curves of the absorber and of the engine to 
which it is to be fixed. The simplified problem, when the torque is applied and 
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the amplitude is measured at /,, will be solved first. The method will then be 
extended to cover the case when torques may be applied at any point, or when the 
severity Of vibration is gauged by the torque amplitude at some other point in the 
system, rather than by the amplitude of displacement of J,. The presence of both 
these effects simultaneously leads to a more difficult problem. 


Writing 2,, for the direct admittance of the crankshaft for torsional motion of 
],, and 2p for the direct admittance of the absorber at the same point, then we obtain 


where 6, is the amplitude of /,. The admittances will be complex quantities, but it 
is assumed throughout that the imaginary part of <,, is negligible; that is to say 
that the damping forces within the engine are sufficiently small to be neglected. 
The maximum values of @, will be those for which 1/2,,+1/ap is a minimum, 
assuming that P, is independent of ». Curves of 1/a¢p and —1/2,, will be plotted 
together, so that, for any given , the distance between the curves will be equal to 
1/2,,+1/ap. The maximum value of @, will then correspond to the value of © 
for which the curves are closest together. 


Because 2p is complex a space curve is required to represent its variation with 
writing 


a 


D 


the values of X, Y, » are then plotted on the rectangular axes as shown in Fig. 3. 
The curve of —1/2,, will lie in the X plane. 


4. Form of the 1/a, Curve 


The damper admittance is given by 


—1 


It is convenient to introduce some auxiliary variables, so that the curve may be 
plotted in dimensionless quantities. Writing 
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which may be written 


K?) , iKk® 
K? 
(K?-1/ + 


This means that for a given ratio of the applied frequency to the natural 
frequency of the damper and for a given value of Q, the quality factor, the reciprocal 
of the admittance may be expressed as a complex multiple of Ap. 


bey 
K 
-xX +X +X 
Fig. 3(a). Fig. 3(6). 
K 


+X 
Fig. 3(c). 


It will be found that X and Y, the real and imaginary parts of 1/ap are 
connected by the relation 


+ =0 


for all values of Q. This shows that the locus of the point (X, Y) in any K plane, 
is a circle of which the centre lies in the »X plane. The circles in all the planes 
form a surface upon which lie the family of curves of 1/2) which is obtained for 
different values of Q. 


If a three-dimensional model of the curves is constructed, it is possible, if 
—1/2,, is first plotted to a suitable scale, to obtain the curves of 1/2,,+1/2p by 
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measurement of the distances between the two curves. Such a model has been made 
in which the circles have been fixed to a transparent base; three wires, which 
represent the curves for Q=2, 5 and 10, have been fixed to the circles. This 
arrangement may be placed on top of a piece of paper upon which the curve of 
—1/2,, for any particular engine has been drawn. The distances between the wires 
and the graph may be measured directly. 


The shape of a typical 1/2» curve is shown by Figs. 3(a), 3(b) and 3(c). 


Fig. 4. Fig. 5. 


Use of the Model to Select a Damper for a System with 
One Degree of Freedom 


The curve of —1/2,, for the system discussed by Timoshenko is shown by the 
full line in Fig. 4. If the natural frequency of the system is fairly close to that of 
the absorber then the existence of two minimum values for 1/2)+1/2,, can be seen 
by inspection of the curves of Fig. 3 and Fig. 4. The selection of a damper which 
will make these two minima as large as possible is the geometrical equivalent of 
Timoshenko’s algebraic method. The ratio /p/J, is first decided because the 
magnitude of Jp is limited by available space at the end of the engine crankshaft. 
For any given ratio of damper frequency to crank frequency a curve of 1/2,, may 
be drawn to the correct scale. It is then necessary to select the damper frequency 
by trial and error, but the procedure may be shortened if a new curve is drawn on 
the base of the model. This curve, shown dotted in Fig. 4, passes through the 
centres of all the circles of the model. It will cut the —1/2,, line at two points 
a and b as seen in Fig. 4. These points give the values of » for which 1/2)+1/2,, 
is independent of »; they are the two points used by Timoshenko. The trial and 
error process now consists of making equal the radii of the two circles of which 
the centres are aand b. It is not necessary to draw a large number of trial curves. 
If one curve is drawn for an engine frequency of, say, 1.1 x absorber frequency, then 
it will give a good approximation to any engine curve corresponding to a frequency 
in the range 1.0 to 1.2 x absorber frequency; it is only necessary:to move the curve 
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bodily. The value of » can then be found by choosing the 1/@p curve for which the 
radii of the two equal circles correspond to minimum values of 1/¢p)+1/2,,. It so 
happens that the two required minima are given by about the same value of yu, 


6. Extension to System with Several Degrees of Freedom 


The —1/2,, curve for the dynamic system of Fig. 1 will be of the form shown 
by the full lines in Fig. 5. Suppose an absorber is to be designed to reduce the 
vibration in the region of the first natural frequency , of Fig. 5. The approximation 
involved in the method discussed consists of replacing the curve of Fig. 5 by 
one of the shape of that in Fig. 4, the simple curve being chosen to cut the vertical 
axis at w, and also to have the correct slope at this point; such a curve is shown 
by the broken line of Fig. 5. It is now necessary to examine how the results of 
this treatment are modified by the difference between the curves of Figs. 4 and 5. 
First it is useful to compare the admittance lines of the anchored one-mass system 
with those of a free two-mass system, Fig 2(b); such a comparison is made in Fig. 6 
in which the line 1 is for the anchored one-mass system. The two-mass 
system lines have an asymptote below the natural frequency »,. If this asymptote 
is much nearer to the origin than it is to »,, then the admittance line in the region of 


3 2 


Fig. 6. 


©, is close to that of the one-mass system; but if the asymptote is closer to w,, then 
the two lines depart from one another appreciably both above and below ,. The 
position of the asymptote is fixed by the ratio of the two moments of inertia; Curve 
2 of Fig. 6 corresponds to a ratio of 5:1 between the disc remote from the absorber 
and the disc to which the absorber is fixed; Curve 3 is for a ratio of 0.5: 1 between 
these two. In the following these discs will be referred to as J, and /, respectively. 
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i A. 


0.8 1.2 1.4 
Fig. 7. 


Now if an absorber tuning is chosen to give two equal circles in conjunction 
with Curve 1, it will be found that with Curve 3 the lower frequency circle will be 
increased in radius and the upper frequency circle will be reduced. This follows by 
consideration of the curves of Figs. 3, 4 and 6. This means that the greatest 
amplitudes above the damper frequency exceed the greatest amplitudes below it. 
This effect will be overcome if the damper frequency is slightly reduced; the 
geometrical treatment shows therefore that if a two-mass system is replaced by a 
one-mass system, and if the Timoshenko design method is then used, the damper 
frequency that is calculated will be too high. 


It has been found by measurement that this effect is not likely to be serious. 
The effect is greatest when the damper mass is small, that is to say when it is about 
10 per cent. or less of the effective engine mass. The effect would also be more 
pronounced if, for any reason, it were necessary to use rather less damping than the 
optimum. In Fig. 7(a)a curve of 2,, 2p/(#,, +%p) against w is plotted for a single-mass 
system. The ratio of absorber mass to system mass is 1: 10 and the tuning, absorber 
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frequency : system frequency, is 1.1:1. This is the optimum tuning as calculated 
by Timoshenko’s method. The damping in the absorber is for Q=5. Fig. 7(b) 
shows the curve which is obtained for a two-mass system, for which the asymptote 
is at 0.8 of the natural frequency, the other data being the same. It will be seen that 
the upper peak is about 20 per cent. higher, but the lower peak vanishes altogether, 
Fig. 7(c) is for the same system as Fig. 7(b) but the absorber tuning is 1.01:1. The 
upper peak is reduced to the first value and the lower peak is present and of equal 
height. Thus there is some advantage in the lower tuning ratio. 


If the system has more than two discs the general conclusions are not greatly 
changed Again, the difference between the admittance curves of the real system and 
those of an equivalent anchored one-mass system depends upon how close the 
asymptotes of the real system are to the natural frequency to which the equivalent 
system corresponds. The presence of an asymptote just below this frequency has 
an exactly similar effect as in the case of the two-disc system. If there is also an 
asymptote just above this frequency then the difference between the two admittance 
curves in this region is reduced. 


7. Variation of Amplitude at Other Points of the System 


The amplitude of vibration at the point of attachment of the damper is not 
necessarily the true criterion of the severity of vibration. For example. in the 
two-mass system which has been discussed previously, the torque amplitude in the 
shaft may be the factor which limits the permissible vibration; the ratio of the 
torque to the amplitude of motion at the damper end will vary with frequency. 
A general expression, in terms of the admittances, may be derived for the amplitude 
at any point of the system when the damper is attached, and the exciting force is 
applied, at another point. Referring to the latter point by the suffix 1 and to the 
former by the suffix 2, the amplitude at 2, due to an applied force P at 1, will be 
Pz,. where z,, is the “cross admittance”; similarly the amplitude at 1 will be Pz,,. 


Both these admittances are to be taken to apply to the system alone, without 
the damper. Now when the damper is attached, the ratio of the amplitude at 2 to 
that at 1 will not be changed. The amplitude at 1 will be obtainable from the 
diagram as before. The amplitude at 2 will be obtained by multiplying by the 
factor 2,,/2,, so that 


Zig 


amplitude = 
P 


Thus if 2,,/2,, is plotted against it will show how the damper frequency should 
be chosen to give the lowest value of the amplitude at 2. 
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ANEW TREATMENT OF THE TUNED ABSORBER WITH DAMPING 


Fig. 8. 


For the two-mass system of Fig. 2(b) 


where the points 1 and 2 are taken to correspond with the discs J, and J,. The 
torque in the shaft will be /,?6, so that if 2,,¢p/(a,,+p) is multiplied by 
Aw? /(A—J, ?) then a curve can be obtained giving the torque in the shaft for unit 
torque applied at 1. The function Aw?/(A—/, w*) is plotted in Fig. 8; the torque 
curve deduced from it, in conjunction with the curve of Fig. 7(c}, is shown in Fig. 9. 
The difference between this curve and that of Fig. 7(c) is very noticeable; the heights 
of the peaks may be made more equal by changing the damper tuning to 1:1, but 
in this case it will be found that the smaller peak is increased and the larger peak 
is not reduced. 


It may be mentioned here that if the magnitude of the exciting force varies with 
frequency, as, for instance, when inertia forces are present, then the amplitude curve 
may be modified in the same manner and again it may be found to be advantageous 
to change the tuning of the absorber. 


8. Excitation at Other Points of the System 


It has been shown how the amplitude at 2 may be found when the exciting 
force and the damper are at 1. The reciprocal theorem, applied to the system and 
damper together, shows that, by the same process, the amplitude at 1 may be found 
when the exciting force acts at 2. Further, if several exciting forces, of the same 
frequency but not necessarily in phase with one another, act at points 2, 3, 
4,.... then the amplitude at 1 may again be found. It will be necessary to find 
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and to obtain thence a curve of the amplitude of 


(S:)p,+ 


where P, is the force at 2, and so on, and where the expression denotes the vector 
sum of the quantities, because P,, P,, and so on, may not be in phase. 


9. More General Problems 


If the absorber is attached at a point 1, the exciting force applied at 2 and the 
amplitude required at 3, then the process becomes more complex. Writing 9, for the 
amplitude at 1, and so on, P, for the applied force at 2, and F for the force exerted by 
the system on the absorber 


6,=P, %3—F 
6,=ap)F 


6,=P, —Fa,,. 


Eliminating F and 6, gives 


3 2 23 


In order to obtain 6, in terms of the measured vector 


1 
{ 


the expression can be rewritten 


6.—P { 2 — + Zo, %3, Zp 
2 
Zi, 


The first term inside the bracket depends upon the system only and is real; the second 
term is a multiple of the distance measured and is complex. 
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A New Approach to Thin Aerofoil Theory 


M. J. LIGHTHILL 
(Professor of Applied Mathematics, The University, Manchester) 


Summary: The general technique’ for rendering approximate solutioas to 
physical problems uniformly valid is here applied to the simplest form of the 
problem of correcting the theory of thin wings near a rounded leading edge. The 
flow investigated is two-dimensional, irrotational and incompressible, and therefore 
the results do not materially add to our already extensive knowledge of this subject, 
but the method, which is here satisfactorily checked against this knowledge, shows 
promise of extension to three-dimensional, and compressible, flow problems. 


The conclusion, in the problem studied here, is that the velocity field obtained 
by a straightforward expansion in powers of the disturbances, up to and including 
either the first or the second power, with coefficients functions of co-ordinates such 
that the leading edge is at the origin and the aerofoil chord is one of the axes, may 
be rendered a valid first approximation near the leading edge, as well as a valid 
first or second approximation away from it, if the whole field is shifted downstream 


parallel to the chord for a distance of half the leading edge radius of curvature p,. 
It follows that the fluid speed on the aerofoil surface, as given on such a straight- 
forward second approximation as a function of distance x along the chord, similarly 
is rendered uniformly valid (see equation (52)) if the part singular like x~' is 
subtracted and the remainder is multiplied by [x/(x+4 p;)}}. 


1. Introduction 


Although the author’s new approach to thin aerofoil theory will be applied 
here only to the two-dimensional flow at negligible Mach number, it is brought 
forward principally in the hope that its use may later be extended to three- 
dimensional flow over wings, and to flow at general subsonic and (in the case of 
swept-back wings) supersonic Mach numbers. In all cases it is only the potential 
flow which is considered, but it is known how this can be used, together with the 
theory of the boundary layer, to predict the real motion of the fluid. 


The existing theory of the two-dimensional incompressible potential flow about 
thin aerofoils (see Goldstein’s Wright Brothers Lecture’ and other references there 
quoted, the basic work being by Theodorsen'*’) is completely satisfactory for the 
practical purpose of obtaining a sufficiently accurate answer in a sufficiently short 


Paper received July 1951. 
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time. However, being based on the theory of conformal mapping, it offers no 
possibility of extension to three-dimensional problems, or to problems involving 
compressibility. 


In contrast, the lifting-surface theory of flow over wings cannot be called 
completely satisfactory, particularly near the leading edge. This theory consistently 
neglects the squares of the disturbances to the main stream, and so is a direct 
extension of the two-dimensional theory called Approximation I by Goldstein (which, 
especially as applied to flow over the camber line only, has a long history, not worth 
recalling here, see, for example, Ref. 1, p. 87). Since in both theories squares of 
disturbances are neglected, the boundary conditions may be satisfied on a plane 
surface, parallel to the main stream, with the same plan form as the wing, because 
the difference between the values of disturbance velocities on such a plane surface 
and on the wing is of the order of the squares of the disturbances. The theory is 
not strictly logical, since for no wing, however thin, is the square of the disturbance 
negligible, at any rate near the leading edge stagnation point, and indeed the 
behaviour predicted in this neighbourhood is manifestly incorrect. 


It is generally believed that the lifting surface theory gives a fair approximation 
to the pressure distribution on a thin wing at moderate incidence, except near the 
leading edge, and part of the evidence for this belief is that the Goldstein higher 
approximations show it to be correct for two-dimensional flow. What is needed is 
a simple method of modifying lifting-surface theory so that it becomes a good 
approximation even near to the leading edge. 


Passing to compressible flow, we may call the theory of two-dimensional 
subsonic potential flow about thin aerofoils, below the shock stall, nearly but not quite 
satisfactory from the practical point of view. The Karmdn-Tsien approach (in its 
practical engineering form, where no deformation of the boundary is contemplated) 
gives results which have been checked satisfactorily with experiment in many cases, 
but its foundations are too obscure to inspire confidence in applying it to any case 
widely different from those checked. The method of Hantzsche and Wendt", on 
the other hand, which is a step-by-step method based on the Prandtl-Glauert theory 
as a first approximation, is very sound and useful where it is applicable, but is 
limited in practice to special shapes, such as Joukowsky aerofoils and elliptic 
cylinders. This is because, although the method involves essentially an expansion 
in powers of the disturbances, it was found necessary at every stage to apply the 
boundary conditions on the surface of the aerofoil, not on a plane surface, owing to 
errors in the higher approximations which otherwise result from the total want of 
accuracy near the leading edge. 


Thus the first approximation, in this method, is not, and must not be, |) 


“ Approximation I” as modified by the Prandtl-Glauert transformation, even though 


this is the solution obtained by consistently neglecting the squares of the | 


disturbances. To get tolerable accuracy in the second approximation it is essential 
to use a first approximation which is valid even near the leading edge. Even in 
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A NEW APPROACH TO THIN AEROFOIL THEORY 


obtaining the second approximation the great simplification which would be 
achieved by applying the boundary conditions on a plane surface is not permissible. 


A parallel difficulty arises in the theory of supersonic flow past wings swept 
back beyond the Mach angle (so that the flow normal to their leading edges is still 
subsonic). In this problem the wave drag cannot be calculated from the ordinary 
completely linearised theory of supersonic flow because of its inaccuracy near the 
leading edge. R. T. Jones) has shown how the wave drag can be corrected, but 
it would be desirable to be able to predict the complete pressure distribution near 
the leading edge, for example in order to detect stalling. 


Thus it is seen that in all these problems involving the flow over thin wings 
with rounded leading edges the mathematically natural method of solution, namely 
expansion in powers of some parameter defining the size of the disturbances, has two 
objections to it: — 

(i) Neither the first nor any subsequent approximation is valid uniformly; 
in fact they are worthless near the leading edge. 


(ii) The inaccuracy due to this increases with successive approximations, so 
that, sufficiently near the leading edge, the process diverges. 


Now in what follows this method of series expansion will be slightly modified, in 
accordance with the author’s general theory of rendering approximate solutions to 
physical problems uniformly valid’, by introducing, with each approximation to the 
velocity field, a slight straining of the co-ordinate system. This is determined 
progressively with the successive approximations to ensure uniform convergence 
near the leading edge. 


Actually, in the problem of two-dimensional incompressible flow which alone is 
considered here, the straining becomes, to a first approximation, a pure translation. 
The velocity field given by the first or second approximation is found to be valid 
uniformly, even near the leading edge, if the field is shifted bodily downstream 
(parallel to the chord) through a distance equal to half the leading edge radius of 
curvature. This radius is itself a multiple of the square of the parameter which 
defines the size of the disturbances, and so the shift makes a difference to the 
equations for the coefficients of successive powers of this parameter in the velocity 
field. As a result of the shift the region of divergence of the process disappears 
altogether, and the singularity, previously at the leading edge, is moved inside the 
aerofoil contour. Here it represents a real feature of the analytical continuation 
of the flow field within the aerofoil contour, of which the dominant singularity is 
approximated by this one at the mid-point of the leading edge radius of curvature. 
But the flow field outside the aerofoil is perfectly regular. To higher approximations 
a real straining, not a pure translation, of the co-ordinate system would be required: 
it might be guessed that this will prove true even to a first approximation when the 
analogous problems of three-dimensional and/or compressible flow are investigated 
on these lines. 
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The work described is novel from the point of view of method insofar as the 
general technique”? is for the first time applied to an elliptic problem. Thus it is 
an isolated singularity which has to be tackled, not a singularity spread along a 
whole characteristic. 

The “natural” process of a series expansion of the velocity field is set out in 
Section 2, and the deduction of a uniform approximation in Section 3. The 
formulae for the velocity distribution on the surface, taking into account the 
incidence, fairing and camber line to a second approximation and uniformly valid 
even near the leading edge, are obtained in Section 4. These are of some general 
aerodynamical interest, even though the problem may be considered previously 
solved for practical purposes. 

For simplicity the author has applied the Kutta-Joukowsky condition at the 
trailing edge, instead of adopting the rather more accurate procedure of using an 
experimentally determined value of the lift-curve slope. This is a modification 
which could easily be made, although it would add to the complication 


of the formulae. 


Notation 


(Dashes Genote differentiation with respect to x) 


Pi __ leading edge radius of curvature 
x distance along aerofoil chord from leading edge 
¢ parameter defining the size of the disturbances 
c  aerofoil chord 
C(x),F (x) functions defining the shape of the aerofoil (see the beginning of 
Section 2) 
y co-ordinate of aerofoil surface perpendicular to chord 
defined by «z=angle of incidence 
u,v velocity components in directions of x and y respectively 
U__ main stream velocity 
a coefficients in the expansion of u in powers of « 
COefficients in the expansion of v in powers of < 
U,,;,U4,;, defined by equations (17) 
U (x),V (x) defined by equations (19) 
W (x),Z(x) defined by equations (22) 
q__ fluid speed at aerofoil surface 
Qmax Maximum value of g 
A,,k,,k, constants in equation (36) 
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he RF real and imaginary parts of 
is nan integer 
a 


Zp»: Optimum value of z, defined in equation (46) 


in p fluid density 

‘he A _ value of u at infinity (in Lemma 1) 

V, value of V (x) atx=0 

li 

al w _ variable of integration 

sly Sarna “s coefficients in the expansions of F (x) and C (x) (see equation (24)) 
the X,Y co-ordinates in the strained system (see equation (26)) 

a X,,X,,..... functions of X and Y in the expansion of x and y in powers of 
ion Vis (see equation (26)) 

ion 


2. Expansion of the Velocity Field in Series 


In this section the problem of determining the two-dimensional incompressible 
potential flow past a thin aerofoil of moderate camber at a moderate angle of 
incidence is treated by expanding the velocity field in a power series in powers of 
a parameter ¢ defining the size of the disturbances. 

To do this let the upper surface of the aerofoil have the equation 
for O0<x<c, and the lower surface the equyation 
y=e[C(x)—F(x)], also for O0<.x<c; where F (0)=C (0)=F (c)=C(c)=0. Then 
the straight segment y=0, 0 << x <c is the chord of the aerofoil as usually defined; 
of the origin and the point (c, 0) are the leading and trailing edges respectively, and the 
equations of the fairing and camber line of the aerofoil are y= +¢F (x) and 
y=C(x) respectively. Now let the angle of incidence be <2. Then the velocity 
components u,v in the x and y directions take the values 


u=Ucosez, v=U sinea ; (1) 


at infinity, while on the aerofoil surface y=<[(C (x)+F (x)] they satisfy (with either 
the upper or the lower sign taken throughout) 


Outside the aerofoil surface they satisfy the equations of irrotationality and 
continuity; in the suffix derivative notation* 


uy—vz=0, ustv,=0. ; (3) 


*uy denotes du/dy, urr denotes (Pu/dx* and so on. 
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If now the velocity field be expanded in powers of ¢ as 


u=U v=U (ev, +2’ 2, (4) 
then by (1) 


u,=0, w,=0, .... v,=0, v, (5) 
at infinity; clearly also the equations (3) are satisfied by each pair (u,,v,) of 
coefficients in (4). But since u,,v, do not depend on « their boundary conditions 
derived from (2) must be expressed not on the surface y=¢ [C (x)+F (x)] itself but 
on a surface independent of <, evidently the two sides y=- +0 of the plane surface 
y=0,0<x<c. Transformation of equation (2) to this surface (which may involve 
analytical continuation of the coefficients u,, v, through the surface of the aerofoil) 
is effected by Taylor’s theorem, giving . 


v+vye(C(xX)tF (xX) +4 +... 
on y= +0, 0<x<c. 


When the series (4) is inserted for u and v on the left hand side of (6) and terms 
of higher order than <* are neglected, it becomes 


{ ev, +2°(v, + F (x)] + + [C (x) + F +4 F } 
x { l-—eu,— (u,+[C(x)+F (x) +22 u,?}. (7) 


Equating this to the right hand side of (6), and isolating the coefficients of ¢, «’, « 
in the resulting equation, we obtain that, on y= +0,0<x<e 


v,=7, (u,t+[C (xX) +F 4, + ry) [C (x) + F (x)] Vay 


a. 


The evaluation of irrotational solenoidal fields (u,,v,) satisfying these boundary 
conditions, with the conditions at infinity (5), depends on a certain well-known 
lemma, which can be interpreted physically in terms of the outflow due to a plane 
distribution of sources. 
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pene I. If u, v satisfy equations (3) outside the strip y=-0, 0<x<c, and if 
u=A, v=0 at infinity, while on the two sides y= +0 of the strip v= +V (x), and 
(4) if u as well as v is integrable from 0 to c for y= +0, then 


u—iv (11) 


(5) X+ly—W 
of It is worth remarking that this lemma, being equivalent to the properties of source 
ns distributions, is essentially trivial, does not depend on the theory of analytic 
Dut functions, and has simple extensions to three-dimensional and compressible flow 
—uce problems. Also it is sufficient by itself to treat all problems arising in the flow at 
Ive zero incidence about a symmetrical aerofoil, which are therefore those for which the 
vil) theory will be most readily extended to other problems (see Section 1). Thus if 
2=0 and C(x)=0, then by (5), (8) and Lemma 1, 
©) u,—iv,= (12) 
zJx+iy—w 
0 
But for dealing with the effect of camber and incidence on u, and 2, it is 
ms necessary to apply the lemma to the function 
X+1y 
—iv,)| ‘ 3 
() where that branch of the complex square root is taken which is positive on y= +0, 
3 0<x<c and negative on y= —0, O<(x<c. This branch is equal to +i at 
infinity. It follows that for a wing with F (x)=0, the function (13) would satisfy 
the conditions of the lemma, by (5) and (8), and therefore 
(8) 
» \t 
= (= { 2+ ) dw}. (14) 
(9) x+iy x+ity—w\c-—w 
For the general case of flow about an aerofoil the complete solution for (u,,7,) 
‘ is evidently the sum of (12) and (14), namely 
“ly 
x+iy-—w x+iy x+iy—w\c-—w 
0 0 
ary | (15) 
wn 
a In this work the position of numerator and denominator in the square root has been 


chosen so that the Kutta-Joukowsky condition of finite velocity at the trailing 
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edge be satisfied. The arguments by which the effect of camber and incidence are 
obtained are not easily extended to three-dimensional problems. The formula (15) 
constitutes Goldstein’s Approximation I for the whole velocity field. 


Letting y—>0 in (15), it follows that 


=U,, (4) +U,, (x) y= +0,0<x<c, . (16) 


- c—w 


0 


where 


and these integrals are to be interpreted as Cauchy principal values. Hence by (7) 
v,=U’ (x)+V’ (x) on y= 4+0,0<x<e,. 
where U (x)=u,, C(x)+4,,F (x) V F (x)4+u,,C (x). (19) 


The problem of finding an irrotational solenoidal field (u,,v.,) satisfying (18) is 
identical with the problem already solved for (u,,v,); the only difference is in the 
conditions at infinity (5). The solution is readily obtained as 


u,—1v,= 2° + + (2%) ) dw 


0 


(20) 


To carry the process one stage further, not because anyone would consider 
using the third order coefficients u,, v, in a practical calculation, but simply for the 
purpose of investigating the analytical nature of the successive approximations, 
(10) is finally written in the form 


v3=W’(x)+Z’ (x) on y= +0,0<x<e, » & 
where 
+ 


0 


+4[C? (x) +F? C” (xX) +C (x) F(X) F” (x), 


. (22) 


0 


+4 [C? (x)+ F? F” (x) +C (x) F (x) C” (x). 
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Then, exactly as before, u, and v, must be given by the equation 


xt+iy—w x+iy x+iy—w \c—w 
(23) 


The nature of the singularities of the coefficients u,, v, at the leading and 
trailing edges will now be considered for n=1, 2, 3. In this investigation it will be 
assumed that the aerofoil contour is an analytic arc except at the trailing edge, 
which is a cusp. This means that C (x) is regular in 0<x<ce (i.e. has derivatives 
of all orders in this interval), and also that x~!(c—x)-# F(x) is regular in O<x<c. 
For example at the rounded leading edge x=0, series expansions of the form 


(F,+F, x+..., x+C,x* +... (24) 


must hold, while at the cusp F must fall to zero at least as rapidly as (c— x), but C 
only like (c— x). 


It is now easy to find the analytic form of u,, v, for n=1, 2, 3 by means of a 
second trivial lemma. 


Lemma 2. If f(w)w!(c—w)} is regular in O<w<c then so is 


i! (w) dw 
x 


(inO0O<x<c). 


0 


It follows from Lemma 2 and (17) that u,, is regular in 0 x<c, and so is 
u,, xt(c—x}t. Hence by (19) both U (x) and V (x)x-!(c—x)-? are regular in 
0<x<c. In this respect U (x) and V (x) are completely parallel to C (x) and F (x), 
and it follows that the singularity of (u,, v,) at either leading edge or trailing edge 
is no worse than that of (u,, v,). Notice how it proved essential to the truth of this 
important first conclusion that C (x) should vanish at both leading and trailing edges; 
otherwise V (x)x~-? (c—x)~? could not have been concluded regular in x<c. 
Thus it is necessary as a first step, before any expansion in powers of the disturbances 
is attempted, that the chord of the aerofoil be selected as x-axis (so that all 


| adjustments of incidence have to be made by changing the velocity at infinity). 


However, the analytic behaviour of (u,,v,) and (u,,v,) is by no means 
satisfactory, at least at the leading edge, where (15) and (20) show that all four 
coefficients become infinite whereas in fact the velocity, even on potential theory, 
must remain finite. Some modification of the theory must be made, and (as Section 
3 will show) it is the analytic behaviour of (u,, v,), now to be investigated, which 


| suggests the modification which is necessary. 
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Now (u;,v3) are given by Equation (23) in terms of functions W (x), Z (x) 
defined in (22). Lemma 2 shows that all the integrals in (21) and (22) are regular 
in 0<x <<, in virtue of the information concerning U (x) and V (x) just obtained, 
It follows that W (x) is regular, using the form of F (x) at the singularities and the 
fact that C(x), as well as being regular, vanishes at both singularities. However, 
Z(x)x-!(c—x)-? is not regular in O<x<c, unlike the corresponding products 
involving F(x) or V(x). The contributions to this product from the first three 
terms of Z (x) in (18), and from the term 4 C?(x)F” (x), are regular. However, the 
term 4 F? (x)F” (x) in Z (x), while of satisfactory order (c— x)! at the trailing edge, 
is of order x-+ at the leading edge. In fact, as a result of this term, 


Z(x)~-}F,°x-? as x—>0. & 


Thus it follows at once by (21) that v, has a higher order of infinity on y=0 as 
x— >0 than v, or v,, namely v,~ + F,* x73. It is clear from this that the 
singularity of the whole vector field (u,, v,) at x=0 is of a similar order, namely 
u, —iv, ~ — py iF (x+iy)-? (although actually the conditions of Lemma | are 
not satisfied and the expression (23) for u,—iv, contains a divergent integral). 


To sum up, there is a single term + 7; F,* x~? on the right hand side of the 
boundary condition for v, ((10) or (21)), in view of which the third order terms 
e*u, and ¢*v, in the series (4) for the velocity field are even more singular at the 
leading edge than the first and second order terms. 


3. Deduction of a Uniformly Valid Approximation 


From the results of Section 2 a uniformly valid approximation will now be 
deduced, following ideas of an earlier paper‘*’, which will first be briefly summarised. 


If there is a singularity near which successive terms in a series solution are each 
still more singular than the last, so that sufficiently near the singularity the series 
diverges, the divergence can be removed by straining the co-ordinate system, and 
hence considering each coefficient of the series solution as a function of new 
independent variables. The required transformation of the independent variables 
may also (like the solution itself) be put in the form of series expansions in powers 
of the small parameter (here ¢), say 


x=X (X, Y)+2? x, (X,Y)+...,° 
y=Y+ey, (X, Y)+2? y,(X, Y)+..., J 


and may be determined progressively as the coefficients in the series for the solutions 


are found. If it is desired only to use the latter series up to and including terms of ; 


order «" then the series (26) may similarly be stopped at this point; here n would 
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be 1 or 2 in practice. The transformation of independent variables alters the 
differential equations or boundary conditions for the coefficients of powers of « in 
the solutions, and can in fact be determined so that the increase in order of magnitude 
of these coefficients near the singularity no longer occurs. In many cases a singularity 
remains but it is pushed outside the domain of physical interest, in which the process 
now converges to a non-singular solution. 


In the problems treated in Ref. 5 the partial differential equation was usually 
hyperbolic, so that the singularity lay along a whole characteristic. The kind of 
transformation required to shift this singular curve (this is necessary, for it is 
essentially a misplacing of the singularity, owing to not letting its position be a 
function of ¢, that causes the divergence) then necessarily took the form of a small 
alteration to that “characteristic variable” which vanished on the said curve. But 
in the present elliptic problem, for which the singularity is confined to a single point 
(the leading edge x= y=0), as it will be generally in elliptic problems, there is more 
freedom in the choice of what variable may conveniently be given a small change 
so as to shift the singularity to the appropriate position, although probably a good 
choice may lead to greater simplicity in the analysis. 


In the present problem it is evident thai a straining of order «? in the 
independent variables will be sufficient since it is only in (u,, v,), not in (u,, v,), that 
an increase in order of magnitude near the singularity first becomes apparent, whence 
there is no need to do anything that will change the field (u,, v,). Further, it is not 
proposed that any more terms in the series for the velocity field than eu, + ¢7u., 
ev, +2?v, be used, and therefore it will be sufficient to terminate the series (26) at 
the terms in <2. Hence all that is required is to determine a transformation 


such that, when u, v are expanded in series of powers of < whose coefficients are 
functions of X and Y, the coefficients of «* are no more singular at the origin than 
those of < and <*. It is assumed without proof, as a result of the careful analysis of 
analogous problems in Ref. 5, that the process could be continued indefinitely and 
that with further, even smaller, modifications to x, y all subsequent increases in the 
order of the singularity could also be prevented. In consequence it is inferred that 
the process converges uniformly within the field of flow. 


The deduction of a suitable form for x, and y, could be made by substituting 
(27) in the equations of motion and boundary conditions, with the coefficients in the 
series (4) for the velocity field regarded as functions of X, Y rather than of x, y, 
and equating coefficients of «, <* and «*. Evidently u,, v,, u, and v, would be 
the same functions of X, Y as they were found to be of x, y in Section 2. But 
u;, ¥, would have a different expression, involving the unknown functions x, (X, Y) 
and y,(X,Y), and these would then be selected so that no singularity in u,, v, of 
worse than inverse square root type occurred at the origin X= Y =0. 
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However, in the present instance it was not necessary to go through this process, 
because the preliminary study of a simple case (the Joukowsky aerofoil), for which 
the exact solution is known, indicated that the simple transformation, obtained 
from (27) by taking y,=0 and x, a mere constant, might suffice. The verification 
that this is true in general was very simple, proceeding as follows. 


Since the co-ordinates are changed by a mere constant shift, the vector fields 
(Un, Vn) must still be irrotational and solenoidal even in the new co-ordinates. The 
conditions at infinity (5) are also unaltered. But the boundary condition (2) takes 
the form 


[= (X, y) 


and when this is expanded in series in the form (6), the left hand side is evidently 
modified only to order <*, while the right hand side is changed, even though terms 
of order <* be neglected, to 


(X)+F’ (X)] + x, [(C” (X)+F” (X)]. 


Hence the equations (8), (9) and (10) remain true, with X written for x 
throughout, except that there is an additional term x, [C” (X)+F” (X)] on the right 
of (10). Hence u,, v,, u., v, are (as foreshadowed above) the same functions of 
X and y as they were previously of x and y. But the boundary condition for », 
in the form (21), again with X for x, is correct only if W (X) be increased by x, C’ (X) 
and Z(X) by x,F’(X). As a result of this change W(X) remains regular in 
0=<X<c, but the part of Z(X) which is infinite at Y¥=0, namely —{F,° X~ 
(see (25) ), disappears provided that 


With the cancellation of this infinite part, what remains satisfies (like F (X) and 
V (X)) the condition that X!(c—X)' Z(X) be regular in 0< X <c, and therefore 
the field (u,,v,) is no more singular as a function of (X,y) than (u,,v,) or (u,,%,). 

Once x, has been determined in this way no further notice need be taken of the 
cumbersome u, and v,. The approximation which has been discovered to be valid 
if the cubes of the disturbances are neglected is 


y)], 
(31) 


where the functions u,,v,, u, and v, were determined in Section 2 by equations (15) 
and (20). A uniformly valid first approximation would still remain if the terms in 
U>, V, in (31) were omitted; but the shift of the origin to the point (4 F,” «?, 0) cannot 
be omitted if the approximation is to have any value at all near the leading edge. 
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This new position of the singularity of the approximate velocity field (31) is 
inside the aerofoil contour, being at the mid-point of the radius of curvature of the 
leading edge. This is evident from the equation of the aerofoil contour 
y=e(+F,x'+C,x+...] and Newton’s formula lim [y*/(2x)] for radius of curvature 
at the origin. 


This justifies the conclusion in the form stated in the Introduction, that the 
velocity field given by the first or second approximation is valid uniformly, even 
near the leading edge, if it is shifted bodily downstream (parallel to the chord) 
through a distance equal to half the leading edge radius of curvature. 


Some departure from the extreme simplicity of this conclusion is to be expected 
when the work is extended to three-dimensional or compressible flow fields. It 
might be found, for example, that a shift in x by an amount depending (at least) on 
z(the distance along the span) was necessary in the three-dimensional case, and that 
it was necessary to take non-zero values for x, and y, in the series (26) in 
compressible flow. Unfortunately the author cannot at present work on these 
suggested extensions, but hopes that some reader may be iriterested enough to 
achieve them, at least in part. 


4. Formulae for the Velocity Distribution on the Surface 


In this section the formulae of Section 3 are specialised to give approximations 
to the distribution of q/U, the ratio of the fluid speed to its value at infinity, over 
the surface of the aerofoil. This is done partly because it makes possible a check 
of the theory against Goldstein’s Approximations II and III, partly because one or 
two new facts of interest emerge, and partly because the simplifications used in the 
process will be helpful as indicating how to proceed when the work is generalised 
as indicated in Section 1. 


On the surface 
q=(u? =u {14+ (C’ F’ } . 


by (2), and either of these forms may be useful in calculating g. It is evident that, 


+ to an approximation in which squares of < are neglected throughout, g may be 


taken as U (1+<u,) and also the surface value of u, may be equated to its value on 
y= +0, and therefore 


q 


on the surface, with u,,, u,. defined by (17). This is Goldstein’s Approximation I. 
Going to the next approximation but without yet making any effort to make it 
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accurate at the leading edge, so that cubes of <* are neglected throughout, we have 


dy +0 [C (x)+ F (x)] +4 ¢? (oy (x) +F’ (x)? 


dw + 


1 | F’(w)+ (w) 
WwW 
0 


0 


on the surface, where use has been made of (15), (20) and the fact that, by (8), 
0u,/dy=0v,/0x=C” (x)+F” (x) on y= +0. It should be noticed that although 
the first and second approximations u,, u,, v,, ¥, to the velocity components on the 
surface are only infinite like x-* at x=0, the magnitude (34) of the velocity on the 
surface, as expanded up to and including terms in ¢’, is infinite like —}«?F,*x~’, 
which arises from the terms <* [F (x) F” (x)+4 { F’ (x) } ?] therein. 


The correction of (33) and (34) near the leading edge is performed by an 
identical process and it will be sufficient here to do it for the more inclusive formula 
(34). The velocity field (15) behaves like 


u,—iv, E (« | (x+iy) 
0 


near the origin, and u,—iv, has a similar behaviour. Thus the velocity field, after 
the translation of axes found to be necessary in Section 3 has been applied, is 
given by 


u—iv 


U 


=1+¢(u, —iv,)+ (u,—iv,) 


=[A,-4i(¢ F,+2?V,)] (x-4F,? +k, F,? 
(36) 


in some neighbourhood of the singularity (} F,? <*, 0) into which the leading edge 
is shifted by the translation; here A,, k,, k,, ... are constants, each containing a 
portion of order « and a portion of order <?. Equation (36) must be approximated, 
for y=¢[C(x)+F (x)], so that it is correct to within O(¢*) for x >>4F,? <*, and 
correct at least to within O(«) for all'x. (The latter condition could hardly be 
improved while the first term of (36) contains the translation of x only to the first 
approximation } F,? 
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e In all the terms in (36) (x—}<«*F,?+iy}! may be replaced by x4+4iyx-}, 
because the squares of these expressions differ only by }(y? x-!—«? F,?)=O (¢? x4) 
which is in the ratio <? to the terms retained when x is not small, and for all x is 
at most in the ratio ¢ to them, this maximum being achieved for x=O(<?). Hence 
u_ gf Ai Vo) — 

+ k,(x+iy)+k,; Gi 


where the higher powers of x!+4 yx! have been replaced by the first two terms of 


*) their binomial expansion, with an error uniformly O (<*). Hence 
), uA, +2 V,) 
oh U t+ x+... + 
he +y{ -4%k, Fk,-3 Fk, xt... 
Now if the cube of « had been neglected without any exception in (38), so that 
the term } y? x~? in the denominator of the fraction had been neglected, the straight- 
forward second approximation 
1+(eu,+¢ so+ye( SH) . ° (39) 
dy y=+0 
35) to u on the surface would have been obtained. Comparing the two (and replacing 
the said } y?x~-? in the denominator by } F,,” «*, which it equals in the only region 
where it is important) we obtain 
—lyx-te 
. (40) 
where in all the terms and (.7 k,)4nx!"~? it has been permissible to multiply 


i by x/(x+4F,? «*) because the error is a factor 1+O(<*) for ordinary x and is only 
dee | 100 per cent. for x—>0, while the terms involved are uniformly O(<). Simplifying 


ga (40) and making use of the value of 0u,/0y on y= +0 as was done in (34), we 
and 
be x4}F22 {(l+eu,+¢é U,)y=+0+ 


first +2*[C(x)+F (x)][C” x-*}. at) 
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It should be noted that this surface value of u is obtained from the “ natural ” second 
approximation by the two simple steps of (ij) removing the term of order x~! by 
adding } F,”<*x~*, and (ii) multiplying by a factor x/(x+4}F,? «?). Thus the zero 
value at the leading edge is achieved exactly. 

Finally, to obtain g/U, (32) shows that into (41) must be inserted a factor, 
which may be written 


to the order of approximation used here, giving 


+ 
4 


Again, this is obtained from the “natural” second approximation (34) by (i) 
removing the term of order x~', this time by adding }¢«*F,?x~-', and (ii) by 
multiplying by a factor [x/(x+4F,?<*)]. These are very simple corrections, and 
it is consoling that the labour involved in obtaining the “natural” second 
approximation would not have to be followed by much additional work to make 
it uniformly valid. 


If the terms in <* were not retained, the approximation (43) would become 


so that equation (33) (Goldstein’s Approximation I) is rendered uniformly valid 
simply by multiplication by a factor. This is in agreement with Goldstein’s 
Approximation II. The factors are not identical, but they agree to within the 
approximation used. 

The calculation of the second approximation (34), which is so easily modified 
to become uniformly valid, would probably be done numerically, in terms of values 
of C (x) and F (x) given at, say 40 points round the aerofoil. The integrals occurring 
in (34), as well as in the definition (17), (19) of U and V, would be computed by the 
convenient Germain-Thwaites-Watson method’. (Of course for a wide range of 
aerofoils given by algebraic formulae, u,,, u,, and hence U, V can be calculated 
analytically, but the part of the integrals in (34) involving U, V would probably be 
too complicated for this approach to be used). 

It is evident from the approximation (44), with the behaviour of (u,)y-+0 at 
x=0 deducible from (35), that near the leading edge the surface velocity is given 
approximately by 


+xt+c! 2 — Sone) (45) 


(x+4F,? 
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fC’ (w) dw 
by where Zopt = wi (c_w)i (46) 
TO 0 
Or, is the angle of incidence for which the leading edge is itself the stagnation point. 
For other angles of incidence the stagnation point is at x=c«?(z—z,,,)’; also the 
peak velocity (if there is one near the leading edge) is achieved at 
12) x=qpy { } , being equal to 
j 
Qmax — U E + (2 | ‘ (47) 
This demonstrates the well-known fact that stalling for ordinary aerofoils sets in, to a 
3) very rough approximation, when the incidence departs from the optimum incidence 

p by the square root of the ratio 4 «* F,*/c of leading edge radius of curvature to chord 
| (this incidence gives Qmax=U ¥ 3). 

y 
nd - The second order formula (34) for surface velocity (which needs modification 
nd as in (43)) is particularly simple if second order terms due only. to the effect of high 
ke incidence are taken into account, for then, by (19), U and V are simply given by 

c—x\t 
ne U (x)= F (x), V (x)= 2( ) C(x). . (48) 
x 
* The resulting formula may be compared with Goldstein’s Approximation III, which 

§ makes no assumption of small incidence, and it agrees to within the order of 
lid approximation. 
ns 
he © Finally the lift of a thin aerofoil is investigated. This is equal to 2xpU times 

| the imaginary part of the coefficient of (x+/y)~* in u—iv at infinity. Clearly this 
” | coefficient will not be affected by the translation of axes. Hence, by (15) and (20), 

the lift is 
1es 
ng 
the 
of 2xpU?| dw |, (49) 
ied 

© §— and the lift coefficient is 

» 
- =u’ (w))(— ) dw . (50) 
en Ce 
0 
Equation (50) gives at once a second order approximation to the no-lift angle, which 
45) 


is not generally known. Also the lift curve slope appears as 


209 


M. J. LIGHTHILL 


c 


w (c—w) 
0 


0 


(51) 


This agrees, to the approximation considered, with the value given by Goldstein’s 
theory when the Kutta-Joukowsky condition is assumed. 


This section, and the paper, may be concluded by giving the second 
approximation to the surface velocity distribution, rendered a valid approximation 
as in (41) even near the leading edge, in a notation free of tiresome epsilons, which 
are really only convenient for the purpose of explaining the theory. If then the 
aerofoil surface is y=C(x)iF (x), O<x<c, where C(0)=C(c)=0 and 
F (x) ~ (2 p, x} as x—>0, so that p, is the leading edge radius of curvature, and if 
@ is the angle of incidence, and if the functions U(x) and V (x) are defined by 
equations (17) and (19), then the ratio g/U of the fluid speed on the surface to its 
value at infinity is given by the formula 


c 


\x+4 pr x-w 


c 


0 


where the upper and lower signs refer to the upper and lower surfaces respectively. 
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The Profile Drag of Yawed Wings of 
Infinite Span* 


A. D. YOUNG, M.A., F.R.Ae.S. 
and 


T. B. BOOTH, B.A., D.C.Ae. 
(Department of Aerodynamics, College of Aeronautics) 


Summary: A method is developed for calculating the profile drag of a yawed 
wing of infinite span, based on the assumption that the form of the spanwise 
distribution of velocity in the boundary layer, whether laminar or turbulent, is 
insensitive to the chordwise pressure distribution. The form is assumed to be the 
same as that accepted for the boundary layer on an unyawed plate with zero external 
pressure gradient. Experimental evidence indicates that these assumptions are 
reasonable in this context. The method is applied to a flat plate and the N.A.C.A. 
64-012 section at zero incidence for a range of Reynolds numbers between 10* and 
10°, angles of yaw up to 45°, and a range of transition point positions. It is shown that 
the drag coefficients of a flat plate varies with yaw as cost A (where A is the angle 
of yaw) if the boundary layer is completely laminar, and it varies as cos * A if the 
boundary layer is completely turbulent. The drag coefficient of the N.A.C.A. 
64-012 section, however, varies closely as cos'A for transition point positions 
between 0 and 0.5 c. Further calculations on wing sections of other shapes and 
thicknesses and more detailed experimental checks of the basic assumptions at higher 
Reynolds numbers are desirable. 


1. Introduction and Review 


The increasing use of sweepback to delay compressibility troubles has aroused 
interest in the behaviour of boundary layers on yawed and swept-back wings, for 
an understanding of this behaviour is required if the characteristics of swept-back 
wings are to be predicted and controlled satisfactorily. This paper is concerned 
only with a limited and relatively simple aspect of the problem, namely the profile 
drag of a yawed wing of infinite span. It is hoped that apart from the intrinsic 
interest of the results, the discussion will add to the basic store of data and ideas 
needed to tackle the more general problem. 


*Part of this paper appeared in the Thesis presented by T. B. Booth for the Diploma of the 
College of Aeronautics, 1949. 


Paper received October 1950. 
[The Aeronautical Quarterly, Vol. IJ], November 1951] 
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The treatment begins with the development of the boundary layer equations 
and the momentum integral equations for laminar and: turbulent flow. The 
development of the laminar boundary layer equations is not new, having been given 
in essentials elsewhere": *:*.“), but it is included here for completeness. It is shown 
that the chordwise flow is independent of the spanwise flow and hence the 
corresponding boundary layer equation can be solved by any of the accepted 
methods. A method for solving the spanwise flow equation is suggested which is 
approximate but sufficiently accurate for the purpose in mind. For a flat plate at 
zero incidence it is found that the spanwise velocity profile is strictly similar at all 
points to that in the chordwise direction, and the resultant direction of flow is 
everywhere that of the free stream. The skin friction coefficient at any point, as 
well as the overall drag coefficient, is reduced by the factor cos! A due to yaw, and 
thus the skin friction at a point is the same as that on an unyawed plate at the 
same distance measured parallel to the undisturbed stream direction downstream 
from the leading edge. This result is peculiar to laminar flow in the boundary layer. 


The turbulent boundary layer equations also lead to the conclusion that the 
chordwise component of flow in the boundary layer can be treated as independent 
of the spanwise component. The momentum integral equation of the former 
component can be solved on accepted lines to yield the component of the drag 
per unit span in the chordwise direction. The drag of unit span of the wing is the 
resultant of this component and the spanwise component. To obtain the latter it is 
assumed that the form of the spanwise velocity distribution in the boundary layer 
is independent of the chordwise pressure distribution and everywhere satisfies a 
power law, of the type normally assumed to hold in the turbulent boundary layer 
on an unyawed flat plate in a zero external pressure gradient. This assumption, 
which might be based somewhat loosely, if intuitively, on the argument that the 
spanwise conditions of flow correspond in some respects to the simplifying conditions 
usually assumed in the theory from which the power law is derived, is shown by 
experiments to be reasonable. For the calculations discussed in this paper the law 
assumed is the +" power law. The spanwise momentum thickness at the trailing 
edge is then derived in terms of the chordwise momentum thickness and readily 
yields the spanwise component of the drag. An expression for the resultant drag 
coefficient is then presented in a simple and concise form. For a fully turbulent 
boundary layer on a flat plate at zero incidence, with the velocity distributions in 
both spanwise and chordwise directions satisfying the + power law, the drag 
coefficient is reduced by yaw by the factor cos* A. Here it may be noted that 
although the flow is everywhere parallel to the undisturbed stream direction the 
skin friction at a point is not the same as that on an unyawed plate at the same 
distance measured parallel to that direction from the leading edge. Were it so, the 
drag coefficient would change with yaw as the factor cost A. 


The method has been applied to calculate the drag coefficient of a flat plate and 
of the N.A.C.A. 64-012 section at zero incidence for various transition positions, 
Reynolds numbers and angles of yaw. In the case of the flat plate (Fig. 2b) the 
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factor describing the variation of the drag coefficient with yaw changes steadily 
from cos* A to cos?A as the transition point moves back from the leading to the 
trailing edge, and there is no noticeable effect of Reynolds number on this factor 
within the range 10° to 10°. In the case of the 12 per cent. thick wing section 
(Fig. 2a) the factor remains very closely cos# A for transition positions from the 
leading edge to 0.5c, and the effect of Reynolds number on this factor is very slight. 


Preliminary experiments made on a flat plate at a fairly low Reynolds number 
(about 10°) tend to confirm the assumptions and conclusions of the theory for angles 
of yaw up to 45°, but further experimental work at higher Reynolds numbers 
is desirable. 


Notation 
The letters L.E. and T.E. refer to leading edge and trailing edge. Suffix “t” 
refers to a time mean. 


x,y curvilinear orthogonal co-ordinates on the surface 
of a body (y is in spanwise direction for yawed | See Section 
infinite span wing) 
z normal distance from the body - 
u,v,w velocity components in the boundary layer in the x, y, z directions, 
respectively 


U,,V, velocity components just outside the boundary layer in x, y 
directions, respectively 


v__ velocity vector=(u, v, w) 
curly 
Q, undisturbed stream velocity 
V 


» components of Q, in direction of chord and span respectively, of 
yawed infinite wing 


A angle of yaw of wing 
u’,v’,w’ turbulent fluctuation components of velocity in x, y, z directions 
respectively 
c chord of wing, measured normal to span 
p ‘pressure 
p density 
coefficient of viscosity 
vy kinematic viscosity =»/p 
D_ drag of unit span of the wing 
X component of D in chordwise direction 
Y component of D in spanwise direction 
Cp[2,8] drag coefficient at a Reynolds number « and an angle of yaw 8 
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Reynolds number=Q, c/v 
boundary layer thickness 


u u 
0 

u 
0 

0 
| Gus 
/Oyy 


component of shear stress in chordwise (x) direction 
component of shear stress in spanwise (y) direction 
value of 7,, at surface 

value of 7,, at surface 

curvature of lines y=constant, x=constant, respectively 
(H,+5)/2 at the trailing edge 

(VV,/U,) G22) 

[—8/(U,)] dp/ ex 

2 (z/5)—2 + 

[z/(66)] 

variable of integration 

[U,/(»x)] 

stream function 

/ (vx 

function defined by v=09/0z 

¢/ / (vx U4) 
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2. Boundary Layer Equations 
2.1. Laminar boundary layer 


The equations of motion of a steady viscous incompressible fluid with negligible 
body forces are, in the vector notation of Ref. 5, 


—VvxX@= — grad (p/p+4v’)—vcurle. (1) 


The equation of continuity is 
G@vv=0. . . (2) 

Consider the flow past a surface and take x, y, z as orthogonal curvilinear 
co-ordinates, such that z is the normal distance from the surface and the lines 
x=constant, y=constant, define parallel networks on the surface. The usual 
assumptions will now be made for the flow inside the boundary layer, namely that 
the rates of change of the velocity components normal to the surface are large 
compared with the corresponding rates of change parallel to the surface and that 
the boundary layer thickness is small compared with a representative linear 
dimension of the surface. It will be assumed further that, if x, and x, are the 
curvatures of the lines y=constant and x=constant, then the quantities «,4, x,8, 
8 Ox,/Ox, 8° Ox,/dy, Ax,/dy, 8? 8° 8° 3° dy, 
and 6° 0*x,/dxdy are all small compared with unity, where 6 is the boundary 
layer thickness. 

With these assumptions the equations of motion yield the boundary layer 
equations : — 


Ou, Ou Ou eu 

ov ov 0% 


The equation of motion in the z direction leads to the result that the pressure may 
be taken as constant across the boundary layer. The equation of continuity is 


Ou dv dw 

dx * dy dz 

For a yawed wing of infinite span, the notation illustrated in Fig. 1 will be 

adopted, i.e. x is measured along a section profile from the leading edge in a plane 

normal to the span and y is measured along the span. In this case the rate of 

change of any quantity with respect to y is zero, and hence equations (3), (4) and (5) 
become 


Ou 1 dp 


dv av 
Aw 
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At the surface u=v=w=0. 
Just outside the boundary layer u=U,, v=V,, du/dz=0v/dz=0, and the flow is 
irrotational, consequently 


(9) 


Since, however, dU, /dy=0 it follows that dV, /dx=0, i.e. V,=constant=V,, where 
V,, is the component of the undisturbed stream velocity along the y axis. 

Equations (6), (8) and (9) and the corresponding boundary conditions are 
identical with those for the boundary layer over the wing in an unyawed flow and 
with a main stream velocity=U,=Q,cosA. Hence, the chordwise flow in the 
boundary layer is independent of the spanwise flow and may be derived by any of 
the well-established methods. An approximate method for solving for the spanwise 
flow will be discussed in Section 3. 


2.2 Turbulent boundary layer 


With the usual boundary layer assumptions the equations of motion for the 
turbulent boundary layer in three dimensions become 


l dp 0 
_ 1ép fay! 4 
p oy oz? ~ Ox (u ( de w (12) 
where undashed letters now refer to time means, the ute denote turbulent 
fluctuations and suffix “t” denotes a time mean. The third equation of motion 
may again be eaeiel as indicating that the variation of pressure across the 
boundary layer can be neglected. The equation of continuity is 


(u’ v’)— (u’w’),, (11) 


(16) 


Again it may be noted from equations (14), (15) and (16) that the chordwise flow 
is independent of the spanwise flow and can be calculated as for the unyawed wing 
in the stream component normal to its span. 
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2.3. Momentum integral equations 


It is usual to assume that the rates of change of the mean turbulent stresses in 
directions parallel to the surface can be neglected, compared with the rates of 
change normal to the surface. Both the laminar and turbulent boundary layer 
equations for general three-dimensional flow can then be written 


(17) 


where 7,,=0u/0Z, 7y,=0v/0z, in laminar flow, and w’),, 
Ov/0z—p(v’ w’),, in turbulent flow. 


Outside the boundary layer 
(19) 


(20) 


If equations (17) and (18) are now integrated with respect to z from z=0 to 
z=6, and use is made of equations (13), (19) and (20), eventually we obtain 


(T22)o= (pU,*6..)+ 8,*)+ Sev, 8,*), (21) 


and ay (p V,? ax (ep V,* Oy2)+ (pV, 6y*)+ (p U, 8,*), (22) 


where (7,.), and (7,.), are the values of -,. and r,., respectively, at the surface, and 


0 0 


is 
(10) ‘ 12 
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For the yawed wing of infinite span, these equations reduce to 


aU, 


0 
(Tz2)o = ax (p U,? 


Writing H,=5,* /6,,, equation (23) can be cast in the familiar form 


[H,+2] 9... 


4. The Calculation of the Profile Drag of a Yawed Wing 


3.1. The drag components 

The drag D of unit span of the wing can be resolved into components X and Y 
along the x and y axis. From the foregoing it follows that X is the drag of unit 
span of the wing when unyawed but with the same transition position in a stream 
of velocity U,=Q,cos A, where Q, is the resultant undisturbed stream velocity 
(see Fig. 1). Hence X can be calculated by any of the established methods (see, for 
example, Refs. 6, 7 and 8). It remains then to determine the spanwise component 
of the drag, Y. To do this, it will be noted that, since the undisturbed stream is 
uniform and the wing section is constant along its infinite span, Y has no form 


Fig. 1. 
Sketch illustrating notation. 
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drag component but is solely the resultant of the frictional stress (7,.),. Thus, per 
unit span of wing, and for each surface 


23) 
Y= AX 
24) 
)dx, f (24) 
= | 5 from (24), 
0 
T.2. 
=p [4,2] 
T.E 
Y where L.E. and T.E. refer to the leading edge and trailing edge, respectively. Here 
nit it has been assumed that the momentum thickness 6,. is zero at the leading edge 
im and is continuous at the transition point, it being accepted that the stress (r,.), 
ity could not become infinite. The problem then reduces to the determination 
or Of [Jee 
nt 
is u v 
now ten 
U 
) 
0 
K6,,, say, 


- (26) 


0 


Suppose for the moment that the value of K is known at the trailing edge of the 
wing. From equations (25) and (26) 


Y=pU,,V, Ky, 


But (see Ref. 6) 


where n= .+5]/2. 
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sin A U, 
Hence Y= SoA Ate. (G* Xx. 


The resultant drag D is given by 


D=XcosA+YsinA, 
and therefore 
D 


Cy [R, A]= =* A+ A sin A 
X { ( U, 
= A A sin? / — 
cos* A +cos A sin? A K, 


\ A { cos? A+sin? A 
[R COS , 0] COS { cos +Ssin* ig) 
(27) 


Here Cy [z, 8] denotes the drag coefficient of the wing at a Reynolds number « and 
an angle of yaw 8. Equation (27) enables the drag coefficient of the yawed wing 
to be calculated in terms of the existing results for the drag coefficient of the 
unyawed wing, the angle of yaw, the chordwise component of velocity just outside 
the boundary layer at the trailing edge and the quantities m and K at the trailing 
edge. The quantity n is known from the chordwise calculations; the problem is 
solved therefore when K,, has been determined. 


3.2 The spanwise flow in the laminar boundary layer 


It is sufficiently accurate to treat the chordwise flow in the boundary layer 
by the Karman-Pohlhausen method or by one of the variants based on it (see, for 
example, Ref. 9). Thus we assume 


+6, (2 +c,(2) +4, (2) 


0? 
and if the boundary conditions » a = ; oe at z=0, 
Ou 
and u=U,,>° = at z=5, 


are satisfied, it is found that -F(2) +G (2) 
8 


-2(3)-2(3) + (3). (8) = 
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8? Op 


and A= 


It then follows that 
37 r 


If it is likewise assumed that 


vol) 


2, 


552 at z=Q, 


and the boundary conditions v 5 


and 


ie. the form of the spanwise velocity distribution is independent of the chordwise 
pressure distributiont. Hence, from (26) 


1 


0 


~ L315 3024) / 315 9072) 
=K (A), sa 

The value of H, is given by 


A A? 


The usual methods of solution for the chordwise flow yield A as a function of x, 
and hence the value of K and H, (and therefore n) can be determined at the trailing 
edge from equations (28) and (29), assuming that the flow remained laminar to the 


tSears® has calculated by a more accurate method the spanwise velocity profiles for a simple 
case of an external chordwise pressure gradient leading to chordwise separation. The spanwise 
profiles show small changes of form up to a point near the separation, after which the change 
of form becomes more marked. These changes are much less severe than the corresponding 
changes of form undergone by the chordwise profiles. 
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trailing edge and was not too near separation anywhere for the foregoing method 
to be seriously in error. In applying equation (27) it is implicitly assumed that the 
solution of the momentum integral equation for the wake developing from completely 
laminar boundary layers follows exactly the lines of the solution for the case when 
the boundary layers are turbulent at the trailing edge. On the face of it this 
assumption seems plausible enough, but in any case the problem of the wing with 
completely laminar boundary layers is at present of little more than academic interest, 


It is noted that for a flat plate at zero incidence with completely laminar 
boundary layers the foregoing leads to v/V,=u/U, or K=1. Since Ur, =U, it 
follows from (27) that 


Cy [R, A]J=cos A Cp [R cos A, 0] 
cos! A Cy [R, 0] (30) 
Cy [R/cos A,0]. 
These results also follow directly from exact solutions of the boundary layer 


equations as can be seen from the following argument. The equations of motion 
are now 


us 
av, 


and the equation of continuity is 


Proceeding along classic lines, from the equation of continuity the existence of a 
stream function y is inferred such that u=0y/0z, w= —0y/0x, and assuming that 
v=f(n)/(vxU,), where »=z/[U,/(vx)] we derive the Blasius equation for 
from the first equation of (31), namely 


ff? +2f’=0. 


Writing v=0¢/0z, and assuming ¢ is of the form ¢=g(n)/(vxU,), the second 
equation of (31) reduces to 


fg” +2g”=0, 
(33) 


with v=U,¢. 
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From (32) and (33) it is seen that 


and hence f’=constant x g’+constant. But f=u/U,, and g’=v/U,, and since 
u=0, v=0, when z=0, and u=U,, v=V,, when z=, we obtain finally 


As already noted, it follows that in this case the resultant flow in the boundary 
layer is everywhere parallel to the undisturbed flow. Equation (30) indicates that 
the drag coefficient is the same as that of an unyawed wing with chord equal to the 
distance from leading edge to trailing edge measured parallel to the direction of the 
undisturbed flow (i.e. c/cos A) in a stream of velocity Q,. This result is peculiar to 
laminar flow in the boundary layer and follows from the fact that for an unyawed 
wing Cp0cl/R*. Any other Cp, R relation would not lead to this result even with 
the flow in the boundary layer everywhere parallel to the undisturbed stream flow. 


3.3. The spanwise flow in the turbulent boundary layer 


To obtain the value of Kr», when the boundary layer is partly or wholly 
turbulent the following assumptions will be made :— 


(i) The distribution of v/V, is the same for all x where the boundary layer is 
turbulent and is independent of the chordwise pressure gradient and therefore of 
the distribution of u/U,, 


(ii) The distribution of 7/V, is then given by v/V,=(z/5)", where m will be 
taken to be +. 


It is highly unlikely that these assumptions describe exactly what happens in 
practice, but experimental evidence tending to indicate that they are sufficiently close 
to the truth for our purposes will be discussed subsequently (Section 5). It may be 
noted, however, that in the usual argument underlying the derivation of the power 
or “log” laws for the velocity distribution in the turbulent boundary layer on a 
flat plate, the development of the boundary layer with distance from the leading 
edge plays no direct part. From this point of view, the argument might be expected 
to apply rather more readily to the spanwise flow in the boundary layer of a yawed 
wing than in the accepted case of the chordwise flow in the boundary layer of an 
unyawed wing. 


With these assumptions we can write 


0 
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From the evidence marshalled in Ref. 7 it may be assumed that u/U, is a function 
of z/6.. and H, only, while @,,/8 is a function of H, only. Hence we can write 


(Us) x, 


6,2 U, 


where K is now a function of H, only, given by 


K=K(H)= | #-[1-(Z) Ja(Z) . 
0 


The function K (H,) can be readily determined using the curves of Ref. 7. The 
value of H, can either be assumed to be a constant over the wing as in the method 
of profile drag calculation described in Ref. 6, or it can be calculated as part of the 
chordwise profile drag calculations as discussed in Ref. 7. Hence the value of H, 
at the trailing edge can be determined, and therefore Ky, can be calculated. 
Sample values of K are as follows: — 


H,=1.29, K=1; H.=14, K=0910; H,=1.5, K=0.777. 


Having determined K and H, (and therefore n=(H,+5)/2) at the trailing edge, 
the drag coefficient can be determined from equation (27). 


For a flat plate at zero incidence it is reasonable to assume u/U,=v/V,, and 
hence K=1. In this case (27) reduces to 


Cp[R, A]=Cp[ReosA, O]cosA. , (35) 
For a fully turbulent boundary layer on the plate the +‘ power law leads to 
Cy[R, 0]=constant/ R* 
Hence Cy[R, O] cost A. (36) 


To illustrate the remarks made at the end of Section 3.2 it may be noted that in this 
case the drag coefficient of an unyawed flat plate with chord=c/cos A would be 


Cy[R/cos A, 0J=Cp[R, cos? A 
ACy[R, A], 


in spite of the fact that the flow is everywhere parallel to the undisturbed 
stream direction. 


4. Specimen Calculations and Results 

The method described has been used to calculate the profile drag coefficient 
of a yawed wing of section N.A.C.A. 64-012 at zero incidences for angles of yaw 
up to 45°, for Reynolds numbers of 10°, 10’ and 10° and for various transition 
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WINGS 
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| |“. | 
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Cp R,0} | 
0:8 | 
FULL LINES REFER TO CALCULATED 4) 
VALUES FOR N.AC.A. 64-012 cos A | \ 
SECTION FOR TRANSITION POINT | | 
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| | 
0:7 
(@) N.A.C.A. 64-012 SECTION AT ZERO INCIDENCE 
A 
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T 
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| 
0-9 
REYNOLDS NUMBER RANGE :- 10° -10° TRANSITION | 
POINT AT 
C,[R,0] | 
0-8 
© EXPERIMENTAL POINTS WITH L.E \ 
WIRE TO PROVOKE TURBULENCE tg 
(R = 1-1 * 108) | 
© 
0:7 


(b) FLAT PLATE 
Fig. 2. 
Calculated values of Cy (R, )/C) (R, 0). 
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positions from the leading edge to 0.5c. Similar calculations, but covering a range 
of transition positions from the leading edge to the trailing edge, have also been 
made for a flat plate at zero incidence. For the aerofoil the calculations were made 
on the assumption that H.. was constant over the section and equal to 1.4. The 
results are shown in Fig. 2 (a) and (b). 


It is of interest to note that although the variation with yaw of the drag 
coefficient of the flat plate is described by the factor cost A when the boundary 
layer is fully turbulent, and changes to the factor cos} A as the transition moves 
back to the trailing edge, the corresponding factor for the wing section remains 
very closely cos! A for transition positions between 0 and 0.5 c. The latter factor 
shows only slight variation with Reynolds numbers. It appears that for a wing, 
where there is usually a small positive pressure at the trailing edge helping to reduce 
the chordwise drag, the factor (U,/U,,)"~' in equation (27) weights the spanwise 
drag component as compared with the chordwise drag component, so that the 
resultant factor Cp[R,A]/Cp[R,0] is somewhat higher than that for a flat plate. 
More extensive calculations are required to throw additional light on the effect of 
section shape and thickness. 


5. Some Experimental Results 


A few preliminary experiments have been made to check some of the 
assumptions of the theory and the results. The experiments were made on a flat 


1-0 
0-8 
uv Az u/U,,ano v/V,, ror = 0 
UV, A= 45° | 
B:- u/U, For @ = sie = 0-885 
C- w/V, ror @ = 32 | 
0:6 
0-4 | 
Ore 0-4 0-6 0-8 10 
7/8 
Fig. 3. 


Measured velocity distributions in chordwise and spanwise planes of boundary layer of yawed 
flat plate at zero incidence and at incidence of 34°. 
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plate of 19 in. chord extending across an open jet tunnel (jet dimensions 
3 ft. 6 in. x 2 ft. 6 in.) of top speed about 125 ft./sec. giving a Reynolds number of 
about 1.2 x 10°. It was intended to check (a) whether the flow in the turbulent 
boundary layer of the yawed plate at zero incidence was everywhere parallel to the 
undisturbed stream direction, (b) whether, with the plate at incidence, the spanwise 
velocity distribution was independent of the chordwise pressure gradient, and (c) 
whether the drag coefficient at zero incidence with a fully turbulent boundary layer 
‘obeyed the law given by equation (35). 


Unfortunately, at the low Reynolds numbers of the tests it was impossible to 
secure a fully turbulent boundary layer or a sharp transition without the aid of a 
turbulence wire of such a size that its drag coefficient was an appreciable and not 
accurately known fraction of that of the plate. The experiments were made with 
the plate at zero incidence and at 34° incidence with various angles of yaw up to 
45°, and consisted of velocity and yawmeter traverses through the boundary layer 
at various stations aft of the leading edge. From these traverses the spanwise and 
chordwise velocity traverses were deduced as well as the drag coefficient of the 
part of the plate ahead of the traverse station. 


The yawmeter traverses made at zero incidence with various angles of yaw 
showed no measurable deviation anywhere in flow direction from that of the 
undisturbed stream, thus checking the assumption that at zero incidence the 
chordwise and spanwise velocity distributions in the boundary layer were similar. 
Some of the spanwise and chordwise velocity distributions measured at x/c=0.855 
at 45° yaw and 0° and 34° incidence are shown in Fig. 3. It will be seen that 
while the chordwise distribution shows a marked change of shape with change of 
incidence, the spanwise distribution shows comparatively little change. The latter 
falls between a 3 and 4" power law for regions in the boundary layer near the 
surface, but fits a somewhat smaller power than 4“ in the outer part of the layer. 


Because of the lack of measurements very near the surface it is difficult to 
estimate the corresponding values of H,=(6,*/4,,) with satisfactory accuracy but 
the value appears to be about 1.4. 


In this connection it may be noted that similar velocity traverses were made in 
the turbulent boundary layer of a yawed wing of finite span and are reported in 
Ref. 10. The wing was of elliptical plan form and of aspect ratio about 5.4: it was 
tested for an angle of yaw of 25° and at angles of incidence of 12° and 14°, at 
which latter incidence part of the wing was stalled. Nevertheless, even under these 
extreme conditions the velocity profiles of the spanwise flow showed themselves 
to be very nearly independent of the chordwise pressure gradients. An indication 
of this is the fact that the value of H, for the spanwise profiles rarely deviated from 
a value within the range 1.3 to 1.5, while the corresponding values of H, ranged 
from about 1.7 to 6.7. The evidence therefore tends to support the assumption 
that the form of the spanwise velocity distribution is relatively insensitive to the 
chordwise pressure gradients. 


227 


i 
ge 
en 
de 
he 
ag 
ry 
es 
ns 
Or 
ce 
se 
he 
te. 
of 

| 
| 

1-0 
ed 


A. D. YOUNG AND T. B. BOOTH 


The measured ratios of the drag coefficients Cp [R, A]/Cp[R,0] are shown in 
Fig. 2(b) for comparison with the theory. It will be seen that the experimental 
results show an even greater reduction of drag coefficient with yaw than is given 
by the theoretical cost A law. This may well have been due to the drag contribution 
of the transition wire, which is almost entirely form drag and would vary as cos’ A, 


It is of interest to note that the theory as developed in this paper leads to the 
following result for thickness at a given chordwise position of the fully turbulent 
boundary layer: — 


A, 


where 6[z, 8] is the boundary layer thickness at the specified position for a Reynolds 
number equal to z and a yaw of angle 8. This relation is shown in Fig. 4 for 
comparison with the measured values. The very close agreement between theory 
and experiment in this instance is probably in part fortuitous, but it is clear that 
the assumptions of the theory cannot be seriously wrong. For comparison the 
factor cos~* A is also shown; this defines what the ratio 6[R,A]/é[R,0] would 


1-4 
CURVE 'A' IS GIVEN BY ASSUMPTION THAT FLOW iS SAME | 
AS THAT OVER UNYAWED WING OF CHORD = C SECA. 


CURVE 'B' IS GIVEN BY ASSUMPTION THAT FLOW IN CHORDWISE } 


DIRECTION IS INDEPENDENT OF FLOW IN SPANWISE DIRECTION. A 
IN BOTH CASES VELOCITY DISTRIBUTION IS ASSUMED TO 
OBEY 1/7th. POWER LAW. 


© EXPERIMENTAL POINTS 


0° 10° 20° 30° 
A 
Fig. 4. 


Effect of yaw on boundary layer thickness on flat plate at zero incidence. 
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be if the development of the boundary layer on the plate depended solely on the 
distance from the leading edge measured in the direction parallel to the undisturbed 
stream and was otherwise unaffected by yaw. 


These experiments cannot be regarded as completely conclusive and further 


work is required, particularly at higher Reynolds numbers, but the results indicate 
that the assumptions made in the theory developed here are unlikely to lead to 
serious errors. 
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Beams on Cross-Girders with Clamped Ends 


J. M. KLITCHIEFF 


(Professor of Applied Mechanics, The Technical University, Belgrade) 


1. Introduction 


In a previous paper''’ it was assumed that the cross-girders are freely supported 
at their ends. In practice these ends can often be considered as clamped (“ built-in”), 
for example, a span of continuous bottom-girders of a ship between two 
transverse bulkheads. 


In the following discussion it is assumed that the ends of the cross-girders 
are clamped, conditions at the ends of the main beams remaining arbitrary, 


as in Ref. 1. 


Notation 


Cc 


m 


total number of lateral loads, or of beams 


distance of cross-section considered measured from left-hand 


end of cross-girder 

length of cross-girder 

length of beam 

moment of inertia of cross-section of beam 

moment of inertia of cross-section of cross-girder 
number corresponding to a particular lateral load or beam 
value of x at point of application of load (see Fig. 1) 
applied load (see Fig. 1) 

total bending moment 

bending moment due to pressure of beam number i 
deflection of cross-girder 


value of u at beam number i 


coefficients of trigonometric series representing deflection 


curve of girder 


Young’s modulus 
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ds n,k,j — integers 
pressure exerted by beam number j on the girder 


Q _ given lateral load acting on one beam 


y.8 coefficients in the relation u,= B(yQ-—£ R)) 


B P/I(ED 

BP Ls(EJ) 
ted Q [m/(8B (B’/B) 
[m/(8 B=) yOR 
[wo 

¢, defined by equation (16) 
lers C, 
ry, 

S C\n+2im| 

and so on 

S Cc”, 
and 


- and so on 


R 
Y 
Vi 
/ 
Fig. 1. 
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R’; _ pressure exerted by beam number i on girders 


land V Assuming 
R”, pressure exerted by beam number i on girders winds sey 
numbered 
Ii and IV 
2R”; pressure exerted by beam number j on central | Ref. | 
girder 


M’; bending moment in first girder due to the pressures R’; exerted 
by the beams i and 12 —i (assuming eleven beams as in Ref. 1) 


=’ summation, in which the terms corresponding to the values of 
n which are multiples of m(n=jm), are not included 


2. Built-in Cross-Girders 


If a girder with clamped ends carries two equal and symmetrically situated 
loads (Fig. 1), the expression for the bending moment 


R| x], for O0<x<&, 


for £<x<L-€, 


r[e(1- £) - for |-E<x<L, 


can be developed in a Fourier series 


— = -(1-cos2nz =) cos 2nx—. 
= n=1 L 


The deflection curve of any symmetrically loaded girder with clamped ends 
can also be represented by a Fourier series 


x 
u 3c, (1 cos 2nz *), ‘ (1) 
n=1 4 


the coefficients C,, being determined by the distribution of the load. Thus, instead 
of equation (14) of Ref. 1, in this case 


( 


and the bending moment due to pressures exerted by a pair of symmetrically situated 
beams, number i and number m-—i will be 


0 


232 


pr 


Ve 


of 


) 


BEAMS ON CROSS-GIRDERS WITH CLAMPED ENDS 


If the total number of beams is an even one, the bending moment due to 
pressures exerted by all of them, is 


}(m- 1) l 1) i x 
M=2% [2 Q2-— 2 (1 cos “Neos 2nt 
i=1 


(1 - cos cos 2n +) cos 2nz * | 
(3) 
Using the known expression for the sum of trigonometric polynomials 
n 4 
l sin 6 
it can be proved that 
s(m—)) j | 0 if mis a multiple of m, | 
>» (1 cos ) (4) 
| 4 m for all other values of n, } 
and that 
0 if m or k is a multiple of m, 
(1 cos cos ) = 1 k=|n+2jm|, Gj=0,1,2...), 
i=1 m m | 
4 m for all other values of k. 
(5) 


Thus, expression (3) becomes 


= x 
DY’ Cinsaim| — cos 4) 
= 


l 
A 2 
<j ~ n 


k 
1 


where the terms corresponding to the values of n which are multiples of m(n=jm), 
are not included in ~’. 


If the total number of beams is an uneven one, the left hand side of (4) must 
be replaced by 


im 
= (1 cos 2nt $(1—cosnrt), ‘ (7) 
i=1 m 
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and the left hand side of (5) by 
( —cos 1 —cos 4 (1 —cos kz)(1—cos nr). (8 
i=1 m m 


It can be proved in a similar manner that these sums are also equal to the values 
given by the right hand side of (4) and (5) respectively. 


Thus the expression (6) holds for any number of beams. Introducing it into 
the differential equation of the deflection curve of the girder it becomes 


where S=3'C,, S,-= =’ Cinteim| » . (10) 
k=1 j 
B. 


With n=1, 2, .... m—1, m+1, m+2, ... equation (9) represents a system of 
an infinite number of linear equations for calculating the coefficients C,. In every 
particular case it can easily be solved by the method of successive approximations. 
It can be solved in the general form also; the corresponding method is developed 
in the Appendix, mainly as a demonstration of the convergency of the series (1). 


3. Several Built-in Cross-Girders. 


The method described in Section 2 will be applied to the example discussed 
in Section 4 of Ref. 1, but assuming now that the ends of the cross-girders 
are clamped. 


Equations (19) of Ref. 1 must be replaced in this case by 


R’;= % 74 2 B (ac 460", 4+ 18C 1 cos kx 


R i= % Q- B 46C’,.+80C k 46C k l coskz 6 (12) 


R 1 —cos kx , 


~ 
k=1 


2 
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and the bending moment in the first girder, due to the pressures R; exerted by the 
beams numbered i and 12 —i, is 


= 26124 1 —cos nt cos 2nz 


44C’,—46C” rt) ( rte | 
( 6C”,.+18C 1 —cos k 6) l—cosn 6 cos 2n L 


Thus, the total bending moment due to the pressures exerted on the first girder by 
all beams becomes 


j 


From the differential equation for the deflection curve of the first girder, 
it follows that 


(13) 


with the notation 


Equation (13) and two similar equations for the second and third girder can 
be written in the form 


6 Bl B’ 1181 
+ S + 46(S” +45”,)+ 1818" +45" N= 
6° Bl 
B il +45) + 80(S” +45») 46 + On 
tpg B +4 4608" + +31 (8 +45 mea 


(14) 


n=1 
of 
S. 
d 

k=1 k=1 k=1 
~ 

d >’ » Ss” a= > n+24) n= C 4 - 
j x / x j x 
rs 
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This system of equations, as system (9), can be solved easily in every particular 
case by the method of successive approximations. For example, for B’/B=0.438, 


416 


OB ” QB OB 


The method for solving equations (14) in the general form is similar to that given 
in the Appendix for system (9). 


APPENDIX 
Solution of the System of Equations (9) in the General Form 
It may be noted that 
and that generally 


Thus, replacing m in equation (9) successively by m+2m, n+4m, ..., |n—2m\, 


|n—4m|, ... and adding all the equations obtained to equation (9), gives 


l + Q) 
where $n 48m* I +cot 2m +3cot 2m (16) 
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Introducing expression (15) into equation (9) it becomes 


Adding together all such equations for n=1, 2, ... gives 


x ~ | 
__! 
n=) (1+42 any (1+42¢,) 


First, the value of S given explicitly by equation (17) must be introduced into 
(15), then the expression obtained for $+ 4S, is introduced into (9), which then gives 
the value of C, explicitly. 


The solution of the problem in its general form can be reduced to the solution 
of three independent systems of equations, the number of unknowns in every system 
being equal to the number of unsymmetrical cross-girders only, instead of being 
equal to the number of intersections, as it would be if the usual method were applied. 
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A Direct Aerodynamic Proof of the 


Biot-Savart Law 


J. LOCKWOOD TAYLOR 
(Institutt for Flyteknikk, Norges Tekniske Hogskole) 


The expression for the induced velocity of a vortex element, 


in Glauert’s notation (where K is the strength of the vortex, h the length of the 
normal, and r the distance from the field point to the element, of length ds) which 
is often known as the Biot-Savart law, and is of fundamental importance in aero- 
dynamics, is frequently left unproved in textbooks. Sometimes reference is made 
to the electro-magnetic analogy, or to a standard treatise such as Lamb. One 
well-known work implies that the law is the solution of an integral equation derived 
from the known induced velocity of an infinite straight-line vortex, but this is 
misleading, since there is an infinite number of solutions, of which only one is 
correct. Another textbook gives a lengthy argument starting from the replacement 
of a line vortex by a sheet of sources and sinks, which is probably not too convincing 
to the average student. Lamb himself gives a long and rather involved proof based 
on Helmholtz and Stokes, with appeals to the Theory of Attractions, with which 
few are likely to be familiar. A simple and direct proof from fluid mechanics seems 
therefore worth while. 


If the law can be established for an elementary square or rectangle, i.e. if it can 
be shown that the total induced velocity of a line vortex of this form is the same 
as would be obtained by summing the elements dq (or their components, since q 
is a vector) then it can be regarded as proved for a closed line vortex of any form, 
plane or not, since it is well known that the area or surface enclosed can be cut up 
into rectangles, and that the sum of these rectangles, regarded as line vortices, is 
equivalent to the original contour enclosing them. It is very easy to prove the law 
for a rectangular line vortex, since a direct solution for the velocity field and 
potential can be written down. 


Consider first the velocity potential represented by 
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in the usual Cartesian co-ordinates, with r?=x*?+y?+z*. Near the x-axis, y=z=0 
nearly and |x|=r very nearly, so that 


as for a two-dimensional vortex, but since r is always positive, the sign of » changes 
with x, so that the direction of the vorticity is reversed at the origin. Similarly 
along the y-axis. Hence altogether the expression for » applies to a cruciform 
vortex, with a reversal in sign at the origin, where the two lines intersect. 


Now consider 


o=tan-' 
zr 


where the suffixes 1, 2, 3, 4 refer to the successive corners of a rectangle. The 
outlying parts of the four cruciform line vortices cancel each other, leaving the 
required solution for a rectangular line vortex. It is a matter of simple algebra to 
verify that the induced velocities according to the expression for # just given are 
in fact in accordance with the Biot-Savart law. The general proof for any vortex 
follows, as previously indicated. If it be argued that the proof applies strictly to the 
whole vortex only, and not to the element, this is quite correct: the law itself applies 
only to the vortex as a whole (although it is stated for the element), since the element 
ds can exist only as part of the whole vortex. 


A REVIEW 


Laws of Free Turbulence H. Reichardt. V.D.JI. Forschungsheft 
No. 414, 2nd impression. Deutscher Ingenteur-Verlag, 
G.M.B.H., Dusseldorf. 1951. 30 pp. Illustrated. 12.50 Dm. 


The first impression of this paper appeared in 1942 (see R.T.P. Trans. 1752). 
This present edition includes all the original data with only slight modifications and 
additions. For example, a new chapter on momentum and heat exchange is included. 


Free turbulence problems consist of those cases of turbulent flow in which wall 
effects are not present. The more important cases are those of free jets and wakes. 
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REVIEW 


The analysis of such problems can be divided into either the phenomenological of 
the statistical approach. In this paper the statistical nature of the turbulent flow is 
not discussed. The former approach, which was originated by Reynolds, owes its 
main impetus to the “ mixing length” concept by Prandtl and the later developments 
due to Prandtl, Taylor and Karman. Their respective theories, namely the 
momentum, vorticity and similarity theories, enabled mean velocity and temperature 
distributions to be predicted with reasonable certainty. It was noticeable however, 
that systematic changes from the predicted values, especially in the first and second 
derivatives of the main mean velocity, did occur and such changes could not be 
dismissed as being due to experimental errors. Recent measurements have shown 
that the distributions and magnitudes of the turbulent velocities which are assumed 
in the “ mixing length” theories are often incorrect and misleading. What is more 
important is the fact stated by Reichardt that all the main results in these “ mixing 
length” theories can be deduced from dimensional considerations alone. 


Reichardt introduces a new theory of turbulence which is based on the results 
obtained from careful measurements of mean velocity and velocity gradient under 
conditions of free turbulence. The main objection to this inductive approach is 
naturally the fact that since the theory is based on distributions of mean velocity the 
results can hardly be expected to predict accurate distributions of turbulent velocity 
and to throw light on the physics of the turbulence producing the phenomenon. 


From the hydrodynamical equations of motion and continuity, and assuming 
that the longitudinal mean velocity distribution can be represented by the error law 
function, Reichardt derives expressions for the mean lateral velocity, the turbulent 
shear stress (or Reynolds stress) and the mean transport of longitudinal 
momentum in the transverse direction. Direct comparison with measured results 
leads to a general theorem which states that “The mean flux of the longitudinal 
momentum which wanders in the transverse direction is equal to the mean 
longitudinal momentum gradient in this direction multiplied by a function, having 
the units of a length, which is independent of this direction.” This function, which 
is analogous to the “ mixing length,” but unlike it, has no direct physical significance. 
This theorem, together with certain predictions from the theory cannot be accepte? 
without further experimental verification. 


In conclusion it should be stated that Reichardt’s theory, although in many 
ways more logical than the early mixing length theories, still has many seriom 
objections. At the best it may be regarded as an empirical theory for the 
determination of the distributions of mean velocity under conditions of 
free turbulence. 
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